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Abstract 

We consider the evolution of quasi-free states describing N fermions in the mean field limit, 
as governed by the nonlinear Hartree equation. In the limit of large N, we study the convergence 
towards the classical Vlasov equation. For a class of regular interaction potentials, we establish 
precise bounds on the rate of convergence. 


1 Introduction and main results 

This work is motivated by the study of the time-evolution of systems of N fermions in the mean field 
regime, characterized by a large number of weak collisions. The many body evolution of N fermions 
is generated by the Hamilton operator 


N N 

Hn = + X'^V{xi -Xj) 

j=l i<j 


( 1 . 1 ) 


acting on 

= {q G : q(x^i,... ,a:^Ar) = a^ipixi,.. .,xn) for all vr G Sn}, 

the subspace of permutation antisymmetric functions in denotes here the sign of the 

permutation vr). Due to the antisymmetry, the kinetic energy in (jl.l|) is typically (for data occupying 
a volume of order one) of the order (for bosons, particles described by permutation symmetric 

wave functions, it is much smaller, of order N). Hence, to obtain a non-trivial competition between 
kinetic and potential energy, we have to choose A = Moreover, the large kinetic energy of the 

particles implies that we can only follow their time evolution for short times, of the order (fhe 

kinetic energy per particle is proportional to the typical velocity of the particles is therefore of 

the order After rescaling time, the evolution of the N fermions is governed by the many body 

Schrodinger equation 


iN^Adti^N^t = 


N 1 ^ 


'^N,t 


( 1 . 2 ) 


j=l i<j 

for G L^(M^'^). It is convenient to rewrite do]) as follows. We introduce the small parameter 


e = Ar-V3 
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and we multiply (11.21) by e^. We obtain 




N N 


i=i 


i<j 


'lpN,t ■ 


(1.3) 


Hence, the mean field scaling for fermionic systems (characterized by the factor in front of the 
potential energy) is naturally linked with a semiclassical scaling, where e = plays the role of 

Planck’s constant. Notice that for particles in d dimensions, similar arguments show that we would 
have to take e = in fact, our analysis applies to general dimensions (with appropriate changes 

on the regularity assumptions); to simplify our presentation we will only discuss the case d = 3. 

From the point of view of physics, we are interested in understanding the evolution of the fermionic 
system resulting from a change of the external fields. In other words, we are interested in the solution 
of (II.3D for initial data describing equilibrium states of trapped systems. It is expected (and in certain 
cases, it is even known) that equilibrium states in the mean-field regime are approximately quasi-free. 

At zero temperature, the relevant quasi-free states are Slater determinants, having the form 


lAsiater (^1) • • • j ) 


1 

7m 


dst (/j (s^i)) l<i j<Af 


where {fj}f=i is an orthonormal system in L^(M^). Slater determinants are completely characterized 
by their one-particle reduced density LOjy, defined as the non-negative trace class operator over L^(]R^) 
with the integral kernel 


UJN{x;y)=N j dX2 . . . dXN'l/jSliiter{x,X2, . ■ ■ ,XN)'tpS\a.ter{y,X2, ■ ■ ■ ,Xn) ■ 

A simple computation shows that 

N 

i=i 

i.e. UN is the orthogonal projection onto the A'-dimensional space spanned by the N orbitals fi,..., fjy 
defining V’Siater (we used here the notation |/)(/| to indicate the orthogonal projection onto / € 
L^(M^)). In the language of probability theory, the one-particle reduced density corresponds to the 
one-particle marginal distribution, obtained by integrating out the degrees of freedom of the other 
{N — 1) particles. Slater determinants have the properties that higher order marginals can all be 
expressed in terms of con via the Wick rule (this is, in fact, the defining property of quasi-free states). 

The many-body evolution of a Slater determinant, as determined by m, is not a Slater determi¬ 
nant. Still, because of the mean-field form of the interaction, we can expect it to remain close, in an 
appropriate sense, to a Slater determinant. Under this assumption, it is easy to find a self-consistent 
equation for the dynamics of the Slater determinant. We obtain the nonlinear Hartree-Fock equation 

iedtUN,t = [-e^A + {V * pt) - Xt,UN,t] ■ (1-4) 

Here pt{x) = N~^cjn{x; x) is the normalized density of particles at x € M^, the exchange operator Xt 
has the integral kernel Xt{x‘,y) = N~^V{x — y)coN,t{x',y), and, as before, e = It is easy to 

check that, if u!N,t=o is an orthogonal projection with rank N, then the same is true for the solution 
coN,t] in other words, the Hartree-Fock evolution of a Slater determinant is again a Slater determinant. 
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In [7], it was shown that indeed, for sufficiently regular interaction potentials, the many body 
Schrodinger evolution of initial Slater determinants can be approximated by the Hartree-Fock evo¬ 
lution, in the sense that the one-particle reduced density associated with the solution of m 
remains close (in the Hilbert-Schmidt and in the trace norm) to the solution of the Hartree-Fock 
equation (ll.4p . Previous results in this direction have been obtained in [9]; convergence towards the 
Hartree-Fock dynamics in other regimes, which do not involve a semiclassical limit, has been also 
established in [SllinillBlIl!. 

At positive temperature, on the other hand, relevant quasi-free states approximating equilibria of 
trapped systems are mixed states, described by a one-particle reduced density lon with trcjjy = and 
0 < wtv < 1 (it follows from the Shale-Stinespring condition, see e. g. [20l Theorem 9.5], that every 
such u}]\f is the one-particle reduced density of a quasi-free state with N particles; Slater determinants 
form a special case, with wat having only the eigenvalues 0 and 1). In the simple case of N fermions 
with one-particle Hamiltonian h = —e^A -|- I4xt and no interaction, equilibrium at temperature T > 0 
is described by the Gibbs state with one-particle reduced density 

where the chemical potential /r £ M has to be chosen so that tiujisf = N. If we turn on a mean-field 
interaction, it is expected that equilibrium states continue to be approximated by quasi-free states 
with one-particle reduced density of the form (II.5h . with the external potential I4xt appropriately 
modified to take into account, in a self-consistent manner, the interaction among the particles (for 
results in this direction see, for example, ffans]). 

In suitable scaling regimes, the state of the system at positive temperature is expected to be well 
approximated by an appropriate mixed quasi-free state. Similarly as in the case of Slater determinants, 
mixed quasi-free states are completely characterized by their one-particle reduced density. All higher 
order correlation functions (i.e. all higher order marginals) can be expressed in terms of watQ. For the 
evolution of mixed quasi-free states, we find the same self-consistent equation (|1.4p derived for Slater 
determinants. We observe here that the properties trwAr = ^ and 0 < ujjy < 1, characterizing the 
reduced one-particle density of mixed quasi-free states, are preserved by the Hartree-Fock equation 
(|1.4I) . In [B], it was shown that, for sufficiently regular potential, the many-body evolution of a mixed 
quasi-free state can be approximated by the self-consistent Hartree-Fock equation (|1.4H (also here, the 
convergence has been established through bounds on the distance between reduced densities). 

To summarize, it follows from the analysis of [TIE] that the many-body evolution of fermionic 
quasi-free states can be approximated by the Hartree-Fock equation (II.4p . This holds true for Slater 
determinants (in this case is an orthogonal projection with rank N) as well as for general mixed 
quasi-free states (satisfying only tr:LJN,t = ^ and the bounds 0 < LJN,t < !)• 

In the mean field regime, the energy contribution associated with the exchange term can be esti¬ 
mated as follows, for bounded potentials V: 


^ j dxdyV{x-y)\u:{x]y)\^ 


< 



II^^AfllHS 


<C1, 


( 1 . 6 ) 


where the full energy is of order N (here we used that the Hilbert-Schmidt norrrH of wtv is bounded 
by A^/^). Because of the smallness of the exchange term, instead of considering the Hartree-Fock 


^In general quasi-free states are characterized by two operators on a one-particle reduced density ujn and a 

pairing density a. Here we restrict our attention to states with a = 0; this is expected to be a very good approximation 
for equilibrium states of fermions in the mean field regime considered here. 

^The Hilbert-Schmidt norm of a compact operator A is defined as ||A||hs = tiA*A. 
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equation (11.41) . we will drop the exchange term and study the fermionic Hartree dynamics, governed 
by the nonlinear equation 

iedtUJN,t = [-e^A + {V * pt),ojN,t] (1-7) 

with pt{x) = N~^LO]\f^t{x',x) (a proof of the fact that the exchange term does not affect the dynamics 
can be found in Appendix A of my 

The Hartree equation (|1.7p still depends on N (recall the choice e = N and the normalization 
trcjAT = N). It is therefore natural to ask what happens to it in the limit N ^ oo. To answer this 
question, we define the Wigner transform of the one-particle reduced density LJN,t by setting 

WnA^, v) = j iON,t + y; ® “ y) . (1.8) 

Hence, WN^t is a function of position and velocity, dehned on the phase-space x M^. It is normalized 
so that 

J Wis[,t{x,v)dxdv = e^tTLON,t = 1 ■ 

The Wigner transform can be inverted, noticing that 

ujN,t{x-,y) = ^ J dvWN,t{ ^ . (1.9) 


Eq. (II.9p is known as the Weyl quantization of the function W^.t- Notice that Hcutv,* ||hs = \/A||ITAr,i||2. 

The Wigner transform ITAr,f can be used to compute expectations in the quasi-free state described 
by ujN,t of observables depending only on the position x or on the momentum —ieV of the particles. 
In fact, for a large class of functions / on 


and 


ii f{x)u:N,t = j dxf{x)ujN,t{x-,x) = N J dvdxf{x)WN,t{ 
tr f{ieV) ujN,t ^ ^ j f{v)WN,tix, v). 


In other words, f dvW]\[,t(x,v) is the density of fermions in position space at point x G M^, while 
f dxWN,t{x,v) is the density of particles with velocity v G M^. Notice, however, that W]y,t is not 
a probability density on the phase-space, because in general it is not positive (this observation is 
related with the Heisenberg principle; position and momentum of the particles cannot be measured 
simultaneously with arbitrary precision). 

From (lEZl), we find an evolution equation for the Wigner transform WN,t'- 


iedtWN,t{x,v) = J dyisdtujN,t ® ^ 

= j dy {-A^+^y/2 + ^x-sy/ 2 ) ^N,t (^ + y; ^ “ y) 

+ j dy ((H * pt)ix -I- ey/2) - {V * pt)ix - ey/2)) ujN,t ® ' 


Using -A^+^y /2 A^_^yi 2 = -2jeV^ • Vy and expanding 

{V * pt){x + ey/2) - {V * pt){x - ey/2) • V(U * pt) + O(e^) 
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we conclude, formally, that 


iedtWN,t{x,v) = j dy'^yUJN,t[x + Y'^x-Y)e 

+ * pt){x)j^^ j dyyujN,t + y) + O(e^) 

= -2iev ■ VxWN^t{x,v) + isV{V * pt){x) ■ V^WN,t{x,v) + O(e^) . 

As a consequence, we expect that, in the limit N ^ oo (and hence e ^ 0; recall that e = 

WN,t approaches a solution Wt of the classical Vlasov equation 

dtWt + 2v • V^lVt = V(V * gt) • V,Wt (1.10) 

with the density gtix) = J Wt(x,v)dv (in contrast with Wjy^t, the limit Wt is a probability density, if 
this is true at time t = 0). The goal of this paper is to study the convergence of the Hartree dynamics 
towards the Vlasov equation (jl.lOh . in the limit V —^ oo. 

This work is not the first one devoted to the derivation of the Vlasov equation (ll.lOp from quantum 
evolution equations. In [ISET], the Vlasov equation is obtained directly from many body quantum 
dynamics, starting from the fundamental A^-fermion Schrodinger equation (the Vlasov equation also 
emerges in the V-boson case, if the mean held limit is combined with a semiclassical limit; see m, 
where the dynamics of factored WKB states is analyzed). In [131 fli] , the authors take the Hartree 
equation (|1.7I) as starting point of their analysis, and they prove convergence (in a weak sense) towards 
the solution of the Vlasov equation (|1.10p . Note that the analysis of [131 E! also applies to singular 
interactions, including a Coulomb potential (the analysis was extended to the Hartree-Fock equation 
in HU). 

In HaimlaoHU, the convergence towards the classical Vlasov dynamics is established in 
an abstract sense, without control on its rate. The problem of determining bounds on the rate of 
convergence is not only of academic interest. When considering applications to real physical systems, 
the number of particles N is large but, of course, hnite. Bounds on the rate of convergence are therefore 
important to decide whether N is large enough for the Vlasov equation to be a good approximation 
of the Hartree and of the full many body Schrodinger dynamics. 

Bounds on the rate of convergence of the Hartree evolution towards the Vlasov equation have been 
first obtained in |3]. In this paper the authors obtain the convergence in the Hilbert-Schmidt with a 
relative rate for sufficiently regular initial data and potentials (they require V G 

and that V G T^(M^, (1 + \p\'^)dp)). For smooth potentials, an expansion of the solution of the Hartree 
equation (|1.7I) in powers of e has been shown in [TU (with no control on the remainder) and in (UEj- 

Our approach here is similar to the one of [3]; we consider the solution of the Hartree equation 
(HZD for initial data lon with sufficiently smooth Wigner transform Wn, and we compare it with the 
Weyl quantization of the solution of the Vlasov equation (jl.lOp . with initial data W]s[. We consider 
regular interaction potentials. In Theorem 12.11 and in Theorem 12.21 we establish bounds on the norm- 
distance of the solution of the Hartree equation u}N,t with initial data ujjy and the Weyl quantization 
ujjsf^t of the solution of the Vlasov equation with initial data Wjsf. For every fixed f G K, the relative 
error is of the order e = in the limit of large N. The dependence on N of these bounds 

is expected to be optimal. This expectation is confirmed by the expansion of (U, where the next 
order corrections are constructed (in fact, if we assumed initial data with smooth Wigner transform 
W]y G VF°°’°°(M^ X M^) and smooth interaction potential V G IF°°’°°(M^), the result of Theorem 12.11 
would follow from Theorem 1.2 in (U). 
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In Theorem 12.11 we get convergence in the trace-norm, for very regular initial data. In Theorem 
E21 we bound the Hilbert-Schmidt norm, under weaker assumptions on the regularity oi Wm- The 
strategy to show Theorem 12.21 is similar to the one of [3]; we regularize the initial data, we compare 
the solutions of the regularized Hartree and Vlasov equations and then we establish stability of both 
equations with respect to the regularization. We can improve the bounds of [3] by using the trace 
norm convergence shown in Theorem 12.11 for the solutions with regularized data. The nonlinearity in 
the Hartree and in the Vlasov equation depends on the convolution of the potential with the density 
of particles in space. Differences among densities can be easily controlled through the trace-norm of 
the corresponding fermionic operators (which are bounded in Theorem [2T]). Estimating them directly 
by means of Hilbert-Schmidt norms, as done in [3], leads instead to a deterioration of the rate of 
convergence. 

Notice that, in Theorems [Q and 12.21 we consider the solution of the Vlasov equation for initial 
data which are not probability densities. The well-posedness of the Vlasov equation for such initial 
data can be obtained adapting the arguments of [8]; in Appendix [A] we sketch the proof. 

If we assume additionally that the sequence of initial data ojm has a limit, in the sense that 
their Wigner transform converge towards a probability density Wq, then we can also establish the 
convergence of the Wigner transform IVv,t of the solution of the Hartree equation towards the solution 
of the Vlasov equation Wt with initial data Wq (in this case, the solution of the Vlasov equation is a 
classical probability density, for all t G K). This is the content of Theorem 12.31 

Our bounds on the norms of the distance between the Wigner transform IVv,* and the solution of 
the Vlasov equation Wt (as well as the bounds for the distance between WN,t and the Weyl quantization 
WN,t of the solution of the Vlasov equation with initial data Wn) hold for sufficiently regular initial 
data. In particular, Theorem l2.2l needs W^q G xM^) (with some additional weights; see Theorem 

l2.4l for the precise assumptions). This condition is justified for initial data describing equilibrium states 
of confined fermionic system at positive temperatures. At zero temperature, on the other hand, the 
system at equilibrium relaxes to its ground state, which can be approximated by a Slater determinant. 
Typically, in this case, the corresponding Wigner transform is not regular. For example, the ground 
state of a system of N free fermions in a periodic box with volume one is a Slater determinant with 
Wigner transform 

Wn{x.,v) = N~~^l{\v\ < ( 1 - 11 ) 

where p = N \s the density of the particles (this system is translation invariant; therefore, particles are 
uniformly distributed in the box). Eq. (|l.lip corresponds to the idea that to construct the free ground 
state, we should fill the N one-particle states with the smallest possible energy (by the antisymmetry of 
fermionic wave functions, there cannot be two particles in the same state). If we switch on an external 
potential and a mean-field interaction, it is believed that the ground state can still be approximated 
by a state with Wigner transform of the form (II.lip : the only difference is that now we have to fill 
low energy states locally, according to an effective particle density /?tf that can be determined by 
minimizing the Thomas-Fermi functional 


£t¥{p) = j dxp^/'^{x) + 


dxVextix)p{x) -h - 


dxdyV{x - y)p{x)p{y) 


among all p G L^nL^/^(M^) with ||/?||i = N. The resulting sequence of Wigner transforms Wiq{x,v) = 
V“^l(|n| < cpl^Yi^)) is not in x M^). So, while Theorem 12.11 Theorem 12.21 and Theorem 12.31 

provide a good description of the fermionic dynamics in the mean field limit at positive temperature, 
they cannot be applied at zero temperature. 
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For such initial data, we do not get norm convergence towards the solution of the Vlasov equation. 
Nevertheless, in Theorem El and Theorem 12.51 we can still prove convergence for the expectation 
of a class of semiclassical observables. Semiclassical observables are functions of the multiplication 
operator x and of the momentum operator —ieV; they detect variations in the spatial distribution 
of the particles on “macroscopic” scales of order one and, at the same time, they are sensitive to 
variations of order in the momentum distribution (corresponding to the “microscopic” length 
scale e). 

Let us stress the fact that, to the best of our knowledge. Theorem 12.41 and Theorem 12.51 are the 
first rigorous results concerning convergence from the Hartree dynamics towards the Vlasov equation 
that can be applied to reasonable approximations of ground states. 

In Section E in the remarks following our main theorems, we provide explicit examples of fermionic 
states, constructed with the help of coherent states, approximating ground states and positive tem¬ 
perature equilibrium states of fermionic systems in the mean-field regime, to which our theorems can 
be applied. 

2 Statement of the results 

In order to state our results in a precise form, we need to introduce some norms for functions on the 
phase space {x,v) € x M^. For s € N, we define the Sobolev norm 

WfWh ='^ [ \'^^f{x,v)\‘^dxdv 

l/3|<* 

where /3 is a multi-index, and can act on both position and momentum variables. For s, o G N, we 
introduce also the weighted norms 

11-^11^5= /(I-ku^)“|V^/(x,u)pdxdu 

We are now ready to state our main theorems. In the first theorem, we assume strong regularity 
of the initial data, and we prove bounds in the trace-norm. 

Theorem 2.1. Let V G IF^’°°(M^). Let ujn be a sequence of reduced densities on with 

trujN = N, 0 < ujn < 1 and with Wigner transform Wn satisfying ||IFAr||j|^5 < C, uniformly in N. 

We denote by iON,t the solution of the Hartree equation 

idtuiN,t = [-e^A + {V * pt),ujN,t] (2.1) 

with pt{x) = N~^u!N^t{x]x) and initial data ujn- 

On the other hand, we denote by WN,t the solution of the Vlasov equation 

dtWN,t + 2v ■ VxWN,t = V(V * pt) • VvWN,t (2.2) 

with pt{x) = f dv Wjv^t(x,v) and with initial data Wna = Wn- Moreover, let uJN,t be the Weyl 
quantization ofW]\f,t, defined as in [LW- 
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Then there exists a constant C > 0 (depending on ||y||^ 2 .oo and on sup^y ||WAr||j:^ 2 , but not on the 
higher Sobolev norms o/Wn) sueh that 
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(2.3) 


k=l 


Remarks. 

1) Recall that we use the normalization trcajv,* = In this sense, (12.311 shows that and u!N,t 

are close, in the limit of large N, since their difference is smaller, by a factor e = than 

their trace norms. 

2) The assumption HVFArll/ri < C on the Wigner transform of the initial data is equivalent to 
suitable commutator estimates for the initial fermionic reduced density ujjsf with the differential 
operator V and the multiplication operator x. We begin by noticing that 



dxdv \VxW]sf{x, u)p 



iV-i||[V,a;w]||fis. 


Similarly, we find ||V^WAr||| = N ^||[ic,ww]|I hs- As for the weights in the definition of the 
i^®-norms of Wn, we notice that 


11(1 + < C^iVi(l + x^- s^A^/^unWIs 


for some ^independent constant C > 0 . 

Proceeding analogously, one can show that the estimate IIVPatIIhI — ^ follows from the bounds 


N ^11(1+ x^ - e^A)“/2[ai, [02, [as, [ 04 , [a5,taAr]]]]]||HS < C 


(2.5) 


uniformly in N and for all choices of ai,..., 05 with either a* = x/e or a* = V. 

Therefore the commutator structure allows to quantify the regularity and decay properties of 
the quantum state Wn- Estimates of commutators [x,ujn] and [eV,a;Ar] already played a key 


role in 0 [7] . 

3) The estimate sup^ I|II0 v|Ih| — ^ equivalently, the bounds (12.51) . are expected to hold true for 


fermionic mixed states, describing systems of N particles in equilibrium at positive temperature, 
in the mean-field regime, [ 6 ]. A reasonable approximation for the reduced density of such a state 
is given by the superposition 



( 2 . 6 ) 


/p,(x) = e-3/2e-*P-/-ff(x-r) 


of the coherent states 


(2.7) 








with a probability density M with 0 < M{r,p) < 1 and 


J dpdr M (r, p) = 1 

In (12.7p . the function g is assumed to vary on the (possibly A^-dependent) scale 6 and to be 
normalized so that ||g '||2 = 1- For simplicity, we shall make the explicit choice 


9{x) 


(27r(52)3/4 


( 2 . 8 ) 


It is simple to check that, with the definition (12.bp . one indeed finds that 0 < wat < 1 and 
trwAT = N. 

The smoothness and decay properties of the Wigner transform IFtv of (|2.6p follow from analogous 
properties of the phase space density M{r,p), i.e. 


WWnWhI < C\\M\\h, . 


(2.9) 


In fact, according to the previous remark, to prove (12.9p it is enough to show (12.5h . To this end, 
we notice that 

[x/e,ujN]{x-,y) = j dpdr M{p,r){-iVp)fpr{x)fpr{y) = J dpdr {iVpM{p,r))fpr{x)fpr{y) 
[V,ujN]{x]y) = J dpdr M{jp,r)Vrfprix)fpr{y) = - j dpdr {VrM{jp,r))fpr{x)fpr{y) . 

( 2 . 10 ) 


More generally, using integration by parts, all commutators oi with x/e and V can be 
written as superpositions of coherent states, weighted by derivatives of the phase space density. 
Therefore, (12.9p follows from 


|(/pr,/pv)| = J dx fpr{x)fp^r'[x) = CiV exp | - ^ - ^{p - p'f^ , (2.11) 


for a constant C > 0, independent of N and 5, and from the bound 

||[ai,[o2,...,[aj,a;Ar]...]||^S < [ dpdp'drdr' M{p,r)\\V^ M{p' ,r')\\{fp^rjp',r')\^ 


< C7iV||V^M||2||||V^M||2 < CN\\M\\ 


( 2 . 12 ) 


jt3 ? 
-^0 


for an appropriate multi-index /3 with \/3\ = j. The effect of the operators (1 + x'^ — e^A) 
appearing in (12.Sp can be controlled using the decay of (12.111) and of the probability density M. 

We conclude that, for any probability density M £ x M^) with 0 < M(r,p) < 1 for all 

r,p £ M^, the sequence of reduced densities (|2.6I) is an example of initial data satisfying the 
assumption of Theorem 12.11 


In our second theorem, we relax partly the regularity assumption on the initial data. To reach 
this goal, we start from (12.3p and we apply an approximation argument. In contrast with Theorem 
12.11 here we only get bounds for the difference — djN,t in the Hilbert-Schmidt norm (the Hilbert- 
Schmidt norm of a reduced density is directly related with the norm of its Wigner transform; there 
is no such simple relation between the trace norm of a reduced density and the L^-norm of its Wigner 
transform). 
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Theorem 2.2. Let V G L^(R^) be sueh that 

j\V{p)\{l + \p\‘^) dp < oo . (2.13) 

Let u!n be a sequence of reduced densities on with trujN = N, 0 ^ ojn < 1 and with Wigner 

transform Wn satisfying ||WAr||j :^2 < C, uniformly in N. 

As in Theorem \2.1\. we denote by uJN,t the solution of the Hartree equation \2.1\) with initial data 
ujn and by u}N,t the Weyl quantization of the solution WN,t of the Vlasov equation 112. with initial 
data Wna = Wn- Then, there exists a constant C > 0 depending only on supjy ||W^Ar||i ^2 and on the 
integral i2.13\) such that 

\\uJN,t - WAT.tllHS < CVNe exp(Cexp(C|i|)). (2.14) 

Instead of comparing the solution u}N,t of the Hartree equation with the Weyl quantization ujN,t 
of the solution of the Vlasov equation WN,t, we can equivalently compare WN,t with the Wigner 
transform WN,t of ujN,t- Eq- (I2.14p implies that 

\\WN,t - WN,t\\2 < C'eexp(C'exp(C'|t|)) (2-15) 

If we assume that the fermionic initial data cuw has a Wigner transform Wjsi (with appropriately 
bounded //|-norm) approaching, in the limit of large N, a probability density Wq on the phase space, 
we can also compare the Wigner transform Wv,* of the solution ww,* of the Hartree equation with 
the solution Wt of the Vlasov equation with initial data Wq. In the next theorem, we show the 
L^-convergence of HOv,* towards Wt- 

Theorem 2.3. Let V G L^(M^) be such that i2.1S\) holds true. Let ujn be a sequence of reduced densi¬ 
ties on L^(M^), with trcuw = N, 0 < ujn < 1 and with Wigner transform Wjy satisfying ||H 7 v||jy 2 < C, 
uniformly in N. 

Furthermore, let Wq be a probability density on x with ||IVo||j:^| < oo and such that 

\\Wn— Wq\\i < Ckn,i, and \\Wn— Wq \\2 < Ckn ^2 (2.16) 

for sequences kw,!) i^n,2 > 0 with knj —)• 0 as V —)• oo for j = 1 , 2. 

Let denote the solution of the Hartree equation f2.1\} with initial data ujn and let WN,t be its 
Wigner transform. On the other hand, let Wt denote the solution of the Vlasov equation \2.2\) . with 
initial data Wq. Then we have 

\\WN,t - Wth < C'eexp(C'exp(C'|t|)) +C{kn,i + KAf, 2 ) exp(C'|t|) (2.17) 

Remarks. 

1) Notice that, if \\Wn — WqUi < ktv,! for a sequence —)• 0, and if HHOvU/fa, ||iyo||j :^2 < C 

1 /2 

uniformly in N, then, automatically, \\Wn — H 0 II 2 < i.e. the second condition in (|2.16p 

1 /2 

follows from the first one, if we take tiN ,2 = However, it is often possible to get a better 

estimate on kn, 2 , improving the bound (I2.17p (for instance, in the example discussed in the next 
remark, we find kn ,2 = i^N,i = 


10 














2) An interesting example of sequence of initial data satisfying all assumptions of Theorem 12.31 can 
be constructed again by means of coherent states. As in (12.61) . consider the fermionic reduced 
densities 

LON{x;y) = j dpdrM{r,p)fprix)Jpj.{y) 

with fpr{x) = — r) and with M a probability density on the phase-space, with 

0 < M(r,p) < 1 and ||M||i = 1 and such that ||M||j |^2 < oo. For simplicity, we choose g as 
in (12.81) to be a Gaussian function, localized on the length scale 5 = 6{N), with 6{N) —)• 0 as 
N —^ oo. 

The Wigner transform of ujn, defined as in (|1.8p . is given by 


WNix,v) = 


(27r)^ 

(27r)3 

23/2 


dy WAT (X -h y; X - y ) 


(2 ^j2 )3/2 I dydrdpM{r,p)e^y 

f , 

/ dr dp M (r, p) e e 




2S'^ 




2 

71 


(27re)3 

where, in the last step, we evaluated the integral over y. We find 

23/2 




< 


(27r)3 

23/2 


dxdvdrdpe i e ^ \M {x + 6r,v + ep/5) — M{x,v)\ 


dxdvdrdpe i e ^ 


f 


X / dX[d\r\\{VxM) {x + X6r,v + Xep/d)\ 


+ tIpI \{'^vM) (x + Xdr,v + Xep/S)\ 
0 


<C<5||V,M||i + C'-||V,M||i 


< C 


5 + 


hJ 


l|M|| 


HI 


and similarly. 


\\Wn-M\\2<C 


<5 + 


hJ 


l|M|| 




To optimize the rate of the convergence Wjsf —>■ M (i.e. to make the sequence of initial data 
as “classical” as possible), we choose 5 = (recall that e = From Theorem 12.31 

we conclude then that the distance between the Wigner transform WN,t of the solution of the 
Hartree equation and the solution Wt of the Vlasov equation with initial data given by the 
probability density Wq = M is bounded by 

||VF;v,t - Wth < Ce'/2exp(Cexp(C'|t|)) 


Although in Theorem 12.21 and in Theorem 12.31 the assumptions on Wat are weaker than in Theo¬ 


rem [2Tl we still need VFat G 


with a norm bounded uniformly in N. As pointed out in the 
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introduction, this assumption is typically satisfied for interesting initial data at positive temperature 
(like the ones constructed in the remarks after Theorem l2.dp . but it is not valid for Slater determinants 
approximating the ground state, which are relevant at zero temperature. 

In the next two theorems, we establish a weaker form of convergence for the solution of the Hartree 
equation towards the solution of the Vlasov equation. We prove convergence after testing against a 
semiclassical observable (whose kernel varies on the length-scale e in the (x — y) direction). The 
advantage of these two results, as compared with Theorems 12.11 and 12.21 is the fact that they require 
much weaker assumptions on the initial data; in particular, they can be applied to reasonable and 
physically interesting approximations of the ground state of confined systems (examples of such states 
are constructed in the remark after Theorem 12.51) . 

Theorem 2.4. Let V € L^(M^) be so that 

J I^(P)I(1 + \p\^)dp < oo (2.18) 

Let LOjy be a sequence of reduced densities on L^(M^), with tiLOjy = N, 0 < lon < 1; such that 

ti \[x,ujisf]\ < CNe, tr |[eV,a;Ar]| < C'.Ve (2.19) 

Denote by Wn £ T^(M^ x M^) the Wigner transform of ojn- We assume that 


WWnWw^o 


E 

\0\<i 


Nix,v)\ < C 


uniformly in N. 

Let the solution of the Hartree equation \2.1\) with initial data ujsf. On the other hand, let 

uiN,t be the Weyl quantization of the solution WN,t of the Vlasov equation Vi.^) with initial data Wn- 
Then there exists a constant C > 0, such that 

|tr I < CNe{l + \p\ + |g|)2e^l*l (2.20) 

for all p,q & t G M. 


Notice that the expectation of the observable appearing in (12.201) can also be expressed in terms 
of Wigner transforms. In fact, for any fermionic operator un, we find 

= j dx {x - eq; x) 

= N [ dxdvWN{x,v)HP'^H‘i'^ = NWN{p,q) 


Hence (j2.20l) can be translated into the bound 

WN,t{p, q) - WnAp^ A) ^ Ce{l + IpI + |g|)^e'^l*l 


where we recall that WN,t is the Wigner transform of the solution ujN,t of the Hartree equation while 
WN,t is the solution of the Vlasov equation with initial data Wn- 

If the sequence Wn has a limit Wq, a probability density on phase-space, then one can also compare 
the Fourier transform of WN,t with the solution Wt of the Vlasov equation with initial data Wq. 


12 










Theorem 2.5. Let V G L^(M^) satisfy Ii2.18\} . Let lon be a sequence of reduced densities on 
with trcjAT = A^, 0 < ojn < 1 and such that 


tr |[x,a;Ar]l < CA^e, tr |[eV,a;Ar]l < CA^e 

Denote by Wn G L^(M^ xM^) the Wigner transform ofiON- We assume that ||W0v||v(/i,i < C uniformly 
in N. 

Furthermore, let Wq G x M^) be a probability density, such that 

\\Wn-W4i<kn 


for a sequence kn with kn ^ 0 as N ^ oo. 

Let iON,t be the solution of the Hartree equation 112.1\) with initial data ujn and let WN,t be the 
Wigner transform of ujN,t- On the other hand, let Wt denote the solution of the Vlasov equation with 
initial data Wq. Then we have 


sup 

p><? 


1 

{l + \p\ + \q\V 


WN,t{p,q) - Wt{p,q) 


< C {e + kn) 


Remark. A physically interesting example of sequence of initial data satisfying the assumptions 
of Theorem 12.51 can be constructed also here with coherent states. Similarly to ()2.6p . we consider the 
sequence of fermionic reduced densities 

ujN{x-,y) = J drdpM{r,p)frpix)frpiy) (2.21) 

with a probability density M G x M^), the coherent states 

fr,p{x) = e-^l\-^^-^l^g{x-r) 

and the Gaussian function g{x) = _ We notice that 


[x,ujN]{x;y) = e J drdp{VpM){r,p)frp{x)frp{y) 

[eV,WAr](x;y) = e j drdp {VrM){r,p)frp{x)frp{y) 

Hence, we obtain 

tr |[x,a;Ar]| < A^e||V^M||i, tr |[eV,a;Ar]| < A^e||VrM||i 

Moreover, it is simple to check that the Wigner transform Wn oi tv n satisfies ||HAr||wi.i < C uniformly 
in N and (similarly to the remark after Theorem 12.,SI) . 


\\WN-M\\i<C{S + e/d)\\M\\wm 

Choosing 6 = we find ||lTjv — M\\i < _ Theorem 12.51 implies therefore that the Wigner 

transform WM,t of the solution of the Hartree equation with the initial data (I2.21D is such that 


D,gGR3 (1 + kl + \P\Y 


WN,t{p,q) - Wt{p,q) 


<CeV2eC|i| 
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for all t G M. Here Wt denotes the solution of the Vlasov equation with the initial data given by the 
probability density Wq = M. Notice that the assumption M G is also compatible with M being 
an approximate characteristic function; this observation is important at zero temperature, to describe 
systems at or close to the ground state. 

The rest of the paper is devoted to the proof of our five main theorems, appearing in Sections 
[3][5j Appendix [B] contains an important lemma on the propagation of regularity for the solution of 
the Vlasov equation (ll.lOh . which is used in Sect. [3] and Sect. [H Appendix O on the other hand, 
contains a bound on the propagation of certain semiclassical commutators, which plays a key role in 
Sect. [Hand in Sect. [5l 

3 Trace norm convergence for regular data 

Here we prove Theorem 12.11 Recall that denotes the solution of the Hartree equation 

iedtUJN,t = [hH{t),ujN,t] 

with the Hartree Hamiltonian 

hnit) = -e^A + {V * pt)ix) 

and the density pt{x) = N~^ujN,t{x;x). We introduce the two-parameter group of unitary transfor¬ 
mations V({t-,s), generated by hH{t). In other words, solves the equation 

iedtU{t;s) = hH{t)U{t; s) (3.1) 

withZ^(s;s) = 1, for all s G ffi. Notice that 0). 

On the other hand, is the Wigner transform of the solution WM,t of the Vlasov equation 
(jl.lOl) . We find that UN,t satisfies 

iedtUJN,t = [~£^A, + At 

where At is the operator with the kernel 

At{x]y) = V(V* pt) • (x - y)^N,t{x;y) 

We conjugate now the difference — ^N,t with the unitary operator U{t]d). Taking the time 
derivative, we find 


iedtU*{t; 0) {iON,t - ; 0) 

= W{t-^){[V * puunA- At)U{t-^) 

= W{t- 0) ([V * {pt - pt),^NA + 0) 


where Bt denotes the operator with the integral kernel 


Btix;y) = 


iV*pt)ix)-iV*pt){y)-^iV*pt) 


x + y 


■ (x-y) 


LON,t{x]y) 


(3.2) 


(3.3) 
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Integration in time gives (since, at time t = 0, ujn,o = WAf,o = ^n) 


{ujN,t -ujN,t)l^{t]0) = — f ll*{t;s) [1/ * {ps - ps),u]N,s]l^{t] s) ds 

Jo 

1 /■* 

H-/ hi*[t] s) BsU{t\ s) ds 

Jo 

Taking the trace norm, we obtain 

tr\ujN,t < - [ tT\[V * {ps-ps),djN,s]\ds + - I trl^^lds. 

^ Jo ^ Jo 

We will estimate the two terms in the right-hand side of ()3.5p separately, and conclude by 
Gronwall’s lemma. 


(3.4) 


Estimate of the first term in (|3.5p . We start by considering the first term on the r.h.s, 
To this end, we observe that 

tr |[I/* (ps - p^),a;Ar,s]| < J dz\ps{z) - psiz)\tT\[V{. - z),djN,s]\ 

< IIps - Psili suptr |[I/(2 - .),WAf,s]|- 

Z 

We start by estimating the last term in the right-hand side of (13.6p . We have 

tr |[I/(- - z),0JN,s]\ = tr 1(1 - e^A)"^(l + x^)~^{1 + x'^){l - e^A)[I/(- - z),lvn,s]\ 

< 11(1 - e^A)-\l + x2)-1||hs 11(1 + x2)(l - e^A)[V{- - z),5iv,.]||HS ■ 


(3.5) 
applying 

of (1331) . 

(3.6) 


(3.7) 


An explicit computation shows that 


||(l-e2A)-Hl + x2)-i||Hs <GViV 

As for the operator D := {1 + x^)(l — e^A)\y{z — .),a;Ar^s]) h has the integral kernel 
T>(x; y) = (1 x^)(l - £^A^){V{x - z) -V{y - z)) v) 

= N{I + x^){l - A^){V{x - z) - V{y - z)) J du 


where we used the definition of ujn,s as the Weyl quantization of the solution Wg of the Vlasov equation, 
with initial data Wq. Taking into account the fact that the Laplacian A^, can act on the potential 
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V{x — z), on the function Wg or on the phase y)/^^ we obtain that 

D{x;y) = N{l + x'^)(y{x - z)-V{y - z)) j Wn,s{^—^,v^ hu 

-Ne^{l + x^){IW){x-z) J Wn,s[^^^,v) hu 

- ^(1 + x2) {V{x -z)- V{y - z)) I {A,Wn,s){^,v) dv 

+ A^(l + x^) (y{x — z) — V{y — z)) J Wn,s dv 

- ^(1 + x2) (Vy)(x -z)- I {ViWn,s){^,v) 

- iN£{l + x'^) {W){x — z) ■ J WN,s(—^^,vj V dv 

- ^(1 + x2) {V{x -z)- V{y - z)) I {V,Wn,s){^^,v) dv . 

7 

= : '^Dj{x]y) 
j=i 


We estimate now the Hilbert-Schmidt norm of the different contributions on the r.h.s. of ()3.8I) . To 
control the term Di, we expand 

Di{x]y) = N{l + x^){V{x - z)-V{y - z)) J TyAr,^ 

= N{l + x^){VVm-{x-y) I 

= iN£{l + x^){VV)iO ■ j (V2fT7v,.)(^^,x)e®®’-(^-^)/"du 
for an appropriate ^ on the segment between x — z and y — z. Using the bound 


1 + x^ < 1 + 2 


x + y V £2 /x - ^ 

2/^2 


and the assumption V G W 


2,00 


we get: 


-DiIIhs < CN'^£^ I dxdy 


/' 


1 + 2 


x + y '^2 — y 

y 


2i 2 


{V2Wn,s] 
2 


x + y 


,v)e^^<^-y)/^dv 


< 


CN£‘^ j dXdr [1 + ^ j{V 2 Wn,s){X, v)e^^'^dv 

CN£‘^ j dXdu(l+ x2)2|V2lUjv,s(X,u)p+ cxy y dXdv\VlWN,si^,v)\‘^ 


< CN£^Wn,s\\hI + CN£^Wn,s\\h^ 


(3.9) 
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Similarly, we control the Hilbert-Schmidt norm of the second term on the r.h.s. of (13.81) : 


||^2||hs < dxdr [l + X^ + eV]^ J WN,s{^,v)e^^"'dv 

< CNe^\\WN,sfH0+CNe^\\WN,sfH2 

Proceeding analogously to bound the Hilbert-Schmidt norm of the other terms on the r.h.s. of (13.8p . 
we conclude that 

PIIhs < C^/N 

Proposition IB. II allows us to control the weighted Sobolev norms of the solution of the Vlasov 
equation by their initial values. We obtain 

PIIhs < [e||Wjv||Hi + 

for a constant C > 0, depending on ||WAr||j|^ 2 . Thus, from (13.7p . we finally find 

tr|[V(.-z),5Ar,.]| <Ce^l*IVe[||W^||^i+e||W^||^|+£2||tT^||^3+e3||vp^||^; 

Therefore, from (13.6p : 

tr|[V * {ps -ps),ujN,s]\ < \\ps -Ps||itr|[v(- - z),ujN,s]\ 

* (3 10) 

+ Ce^^^^Ns^Wps - PsWiW [WWnWhI + e\\WN\\Hi + 

= I + II 


Consider first I. We have 
ll/Os - /Oslll = sup 

JgLoo(R3),||J||^<l 

where on the r.h.s. the supremum is taken over all bounded operator with operator norm lesser or 
equal than one. We conclude that 


J{z){psiz) - Ps{z))dz 


< N 


-1 


sup |tr J(cJAr,5 - CJAr^s)! 
JilIJIKl 


IIPs - /Oslll < 4r IwAT^s - 


Therefore, 

I < Ce^l^'etrlo;^,, -w,v,.|||kT7v||Hi (3.11) 

To bound II, we write: 

IIPs - Pslli < \\Ps\\i + IIPsIIi = N~'^tiu;N,s + ll/Oslli < 1 + llHW.slli 
Using that the Vlasov dynamics preserves the U’ norms, we get: 

||IU7V,s||l = IllVivlll = 11(1 + + u2)-2(l + + u2)2W,v||l < C\\Wn\\hO 
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and thus: 


(3.12) 


II < Ce^\^\Ne^\\WN\\Hl [\\Wn\\hI+4Wn\\hI+^“^WnWhi] 
From (13.60 . (I3.10p . (13.lip . (I3.12p . we obtain: 

■y rt rt 


£ Jo 


tr \[V * {ps - Ps),u}N,s]\ds < C I |wAr,s - 

Jo 




\\WN\\H^+e\\WN\\H^+e^\\WN\\m 


(3.13) 


where the constant C > 0 depends on IIVF/vIIh^) but not on the higher Sobolev norms of IFat. This 
concludes the estimate of the first term in the right-hand side of (|3.5p . 

Estimate for the second term in (13.5|) . To conclude and apply Gronwall’s lemma, we need to 
bound the second term in (13.5p . We find 

tr|5,| < ||(l-e2A)-Hl + x2)-i||Hs ||(l + x2)(l-e2A)5,||Hs < CVN\\il + x^){l-e^A)B,\\ns (3.14) 

Let Us := V * ps- The kernel of the operator B := (1 — e^A)Bs is given by 


B{x-y) = N Us{x)-Us{y)-VUs 


X + y 


-Ne^ 


Ne^ 


1 


AUs{x) - 7 AVI/, 


x + y 


{x-y) 

1 


{x-y) - -AUs 


Wn. 
x + y 


V e 


Us{x)-Us{y)-VUs 


x + y 


+ N 


Us{x)-Us{y)-VUs 


^VUsix) - 


2 

-Ne 

-Ne 

7 


1 . 


VUs{x) - -V^U 


x + y 
2 

x + y 
2 

x + y 


■{x-y) J (AiITAr,s) 

{x-y) [ 


x + y 
2 ’ 

I Wn,s 

x + y 


(x — v) 


x + y 


-,v e 


(x — v) 

^^edv 


(x — v) 

-^^edv 


V e 


rx • (x — v) 

,V\V^U^-edv 


{x-y)- VUs 


{x-y)- VUs 


x + y 
2 

x + y 


(ViITjv,s) 

x + y 


Wn. 


x + y 
2 ’ 

,v)ve 


V e 


(x — v) 

^^■ edv 


(x — v) 

^^■ edv 


Us{x)-Us{y)-VUs 


x + y 


{x-y) {v-ViWn,s) 


x + y 


-,v e 


(x — v) 


J2Bj{x;y) 

i=i 


(3.15) 


In the contributions Bi, B 4 , Bq, Bj, we need to extract additional factors of e; the goal is to show that 
11(1 -|- x^)i?||HS < C'/Ne^. To this end, we write 


Us{x)-Us{y)-VUs 


x + y 


{x - y) 


= ['dX[VUs{Xx +{I-X)y)-VUs{{x + y)/2)] - {x-y) 

Jo 

3 el 

J dX J dpdidjUs(^p{Xx + {1 - X)y) + {1 - p){x + y)/2^{x - y)i{x - y)j(^X 


* ,1=1' 
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and we estimate, using the assumption V G and integrating by parts, 


~ ^ r 

\Bi{x-,y)\ <CNe‘^ J 




^ ^ y y\^iv-{x-y)/e 


*,i=i 

Hence, proceeding similarly as we did in (13.9p . we get: 

3 


{l + x^)B,rm<CNh^ 


*J = 1 


dvAjWN,s 




-,v}e 


= CNe^ J dXdr [l + + e \^]" J d,,dy. Wn,s {X, v) e'"" 

< CNe^\\WN,s\\l2 + CNs^\\WN,sfm 

The Hilbert-Schmidt norm of the other terms in the right-hand side of (I3.15p can be estimated in a 
similar way. To do this, it is useful to notice that: 

llV^HlIoo = II w * Vp.lloo < II W||oo||Vp,||i < 

where we used that V G and that ||V;Os||i < C'||lTAr^s||j:^i. The final result is: 

11(1 + X^)i?||HS ^ C^/N e^||HAr^s||j:^2 -|- ||1T7 v,s|| j|f3 -I- e^||lTAr^s||_H'| + S^II^A^jsIIhI 

Therefore, by Proposition IB. 11 

11(1 +a:2)H||HS < [||1 T;v||h| + e|l^iv||/^3 + 

where the constant C > 0 depends on ||HAf||H| not on the higher Sobolev norms of IT/v. This 
gives: 


triHJ < Ce^l'liVe^ 


IIITtvIIj^? +e||lP^||„3 +£2||iy^|| 4 +e3||lT^||^5 


(3.16) 


Proof of Theorem 12.11 We are now in the position to conclude the proof. Inserting (13.131) . (|3.16l) 
into P3.5I) . we get: 


tr \uN^t - ^N,t\ <C tr |wAr,s - a;Ar,s|ds 

Jo 


+ Ce^l*liVe 


l|W^ivllH2 +esup||lTAr||„3 sup ||lTAr||^4 -b sup ||lTAr|| 

^ N ^ N N 


Finally, Gronwall’s lemma implies the desired bound 
tr |wAr,t - ujN,t\ 


< C'A^eexp(C'exp(C'|t|)) 
with C depending only on HITtv 11 /^2. This concludes the proof. 


SUp||lTAr||„2 -besup||lTAr||„3 -be sup ||lTAr|| -be sup ||lTAr|| 
, N * N ^ N ^ N 


□ 
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4 Hilbert-Schmidt norm convergence 


Here we prove Theorem 12.21 and Theorem [231 The proof of Theorem l2.2l is based on an approximation 
argument, together with our previous result Theorem 12.11 


Regularization of the initial data. We start by approximating the initial data Wm- For A; > 0, 
we define 


gk{x,v) = 


(27r)' 




and 


W^{x,v) = {Wn * gk){x,v) = j dx'dv'gk{x - x',v - v')Wn{x',v') 
Then, we have ||lT^||j|^5 < oo for all V G N. In fact, we find 


\W^\\Hi<C\\WN\\H2 ifi<2 and 

for j = 3,4,5. 


Furthermore, we notice that 




(4.1) 


(4.2) 


for s = 0,1 (with the convention = L^) and for a < 4. We denote by the Weyl quantization of 


W^. We observe that 


uj%{x-,y) = N J dvW^ 




J dvdx'dv'e 2 ( 2 “ *) g ^^lTAr(x', n')e*''" 


Nk^ 


(27r)3 

t/ 


(27r 


, x-y , x-y\ 

\ X ^ -- Y~ ^ 

x — y 


(4.3) 


(27r)3 / 


dwdz e yx + ^ 

dwdz e-^'/2g-«;V2 


X z 


^.v —2=-V -iw-- 


VkeeVk ' ojNe ^ e Vke 


{x]y) 


Hence ui^, as a convex combination of fermionic reduced densities, is again a fermionic reduced density 
(i.e. 0 < uj^ < 1 and trcj^ = N). From (14.2p . we hnd 


||a;iv - w^IIhs = ^/N\\Wn - IF^Ib < ^J^\\WN\\H^ (4.4) 

We denote by ujN,t and the solution of the Hartree equation with initial data un and, respec¬ 
tively, uj^. On the other hand, djN,t and will denote the Wigner transform of the solutions WN,t 
and of the Vlasov equation with initial data Wn and, respectively, IT^. Notice that, since the 
Vlasov equation preserves all the LP norms, ||ui7v,t||HS = -V^/^||IFAr^ 4 II 2 = V^/^||lTAr ||2 and, similarly, 
P^,JIhs = II 2 , for alH G M. 
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We need to compare with To this end, we will first compare with Later, we 

will have to compare with and, separately, with 

Comparison of with To begin, we prove that there exists a constant C > 0 such that 

3 

\\^N,t - ^N,t\\iiS < CN^/'^e exp(C'exp(C'|t|))[l + (4.5) 

i=i 

The constant depends on snp^ ||VLv||j|^ 2 , but not on the higher Sobolev norms. To show (14.51) . we 
shall use our previous result. Theorem 12.11 In fact, from (|2.,3p . (14.ip we find 


\\^N,t - ^N,t\\tr < CNe exp(Cexp(C|t|)) \\W^\\hI + X] Wn\\hP+^ 

\ ^ ; TV 4 




13=1 


(4.6) 


< CNe exp(Cexp(C'|t|)) j 1 + '^{eVky 

i=i 


for a constant C > 0 depending only on sup^y ||IL 7 v||jy 2 . We shall use this result to prove an estimate 
for the Hilbert-Schmidt norm of the difference of the two evolutions. Proceeding as in (13.ip - (13.5p . 
we have: 


II‘^T- 5 T|Ihs<) r* [F.(pJ-®, 5 S,,.] 

^ Jo 


HS 


+ 


ds ||i?s IIhs 


(4.7) 


where Bg is the operator with the integral kernel 


Bs{x-,y) = 


V - V *rs{y) - y{y *rs) 


x + y 


{x - y) 


^nAx^v) 


We shall estimate the two terms in (|4.7I) separately. We start with the first. We have: 

[r.(p;'-S).2Tl||„5< / <iz|4W-pi(2)l-|||r(2-.),2Tl 


<llrf-aili/<ip|r(p)| 


HS 


Using \\p^ - p^lli < N the identity 


Jo 


ip-X,U]j^ 




and the assumption (I2.13p on the potential, we conclude that 


[V*{p':-rs)Aks] < CN-y\u;As-djks\\A\[xAks]\\m 

Hb 


(4.8) 


We shall use the regularity of to extract a factor e from the commutator in ()4.8I) . We have: 


[x,^NAA\y) = {x-y) dvW^ 


k (x + y _^\jy.3ezy 


V e 


= e (ivVyW^g 


k (X + y \av^ 


V e 
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and thus, similarly to (|2.4p . 

The second inequality follows from the propagation of regularity for solutions of the Vlasov equation, 
proven in Proposition IB. 11 Inserting the last bound and (I4.6p in ()4.8p . we obtain 


* {ps - Ps),^N,s] exp(Cexp(C'|t|)) j 1 + ^(e\/fe)J 

i=i 


(4.9) 


which concludes the estimate for the first term in (14.7p . Let us now consider the second term in (j4.7l) . 
We have: 


I-^sIIhs = / 


/' 




X + 2/ 


{x - y) 


|wiro(a;;y)P 


J dxdy\uj%^^{x-y)\‘^\x - 


dX 


V(V * p^)(Ax + (1 - X)y) - V(V * ^)((x + y)/2) 
Jo 

< j dxdy\uj%^,{x-y)\‘^\x - yl"^ 

X [' dX r dfi (A - 1/2)v 2(V * pJ)(/i(Ax + (1 - X)y) + (1 - ;r)(x + y)/2) 

Jo Jo 

<C J dxdy\x - y\'^\u}%^^{x-y)\‘^ 

using the assumption (|2.13l) . Since 


7' 




[x - y) y) = -£ I dv 

we find, similarly to (12.41) . 

\\Bs\\hs < CNe^WAyW^jl < C'e^l*le^V||Wi^||^2 < 

where we used again Proposition IB.11 This concludes the estimate of the second term in (HU). 
Therefore, plugging the estimates (14.91) . (|2.4p into (|4.7p . we have; 


\\^N,t - ^N,t\\ns < exp(Cexp(C'|t|)) (1 + ^(e\/fc)^) 

i=i 


as claimed. 

Comparison of ^ with The next step is to compare the Hartree dynamics of the regularized 
initial data with the Hartree dynamics of the original data. Our goal is to show that: 

(4-10) 

for a suitable constant C > 0, dependent on supjy ||hPv||_ff| but not on the higher Sobolev norms. 
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Let U{t\ s) be the unitary group generated by huit) = — + V * pt, with pt{x) = N x). 

From 

iedtU*{t-, 0 )uj%^tU{t- 0 ) = -U*it- 0 ) \v * {pt - p’:),oj%2 ^(i;0) 


we have: 


^N,t - = U{t; 0) (wjv - U*{t, 0)u!%^tU{t-, 0)) 0) 


1 


= U{t-, 0 ){ujN - 0J%)U*{t-,0) + - / dsU{t;s) V * {p, - p';),u;% , U*{t-s) 

le Jq I ’ J 


Hence 


||wAr,t < II^AT -W^IIhS + ^ 


(4.11) 


We start by estimating the commutator in the right-hand side. We have 

Jo 

By Proposition lC.il it follows that: 


ll[e‘’"“’.‘^ClllHS < IpIIIIi.i^CIIIhs < Cble^l*' (l|[x,a.S,]||Hs + ll[eV,i..^]||Hs 


Since 

||[x,a;^]||HS = eN^/^\\V,W^h < £N^^^\\Wo\\h^ 

||[eV,u;i,]||HS = eiV'/'||V,Wo"||2 < eN^/^WWoWh^ 

we conclude that 

ll[e*'’'",wL]||HS<CAV2e|p|eCld (4.12) 

Then, we are left with estimating the trace on the right-hand side of (|4.11l) . To do this, we shall use 
the following lemma. 


Lemma 4.1. Under the same assumptions of Theorem IH.HL there exists a constant C > 0, only 
depending on sup^y but not on the higher Sobolev norms, such that 


sup 


pma 1 + IpI 


tTUP-^{u;N,t - u;%^t)\ < Ce^^'^N 


(4.13) 


Plugging (14.121) . (14.131) into (14.lip , and using the bound (|4.4I) on the difference of the initial data, 
we get 

\\^N,t - w^,tl|HS < 

which concludes the proof of (I4.inp . Thus, we are left with the proof of Lemma l4.11 
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Proof of Lemma \4.1\ Consider, for an arbitrary p G 

tr — ^N,t) — 0)e^^'^U{t] 0 ){ujn — U*{t, 0) ui^ ^U{t] 0)) 

where, as in (1,3.ip . U{t]0) denotes the unitary group generated by hnit) = —e^A + V * pt, with 
Pt{x) = x). From 


iedtU* {t-0) u;% = -W{f,0) V * {pt - pi), U{f,0) 


we find 


= tT:W{t]0)e^P'^U{t-,0){ujN - 0J%) 

-- [\rWit-s)e^P--Uit;s) \v * (p^ - p1),cv%^, 

Jo L 

= trU*{t; 0)P‘^'^U{t-, 0)(a;Ar — 


ds 


+ 


hi! J {psiP)-f^siP)) trU*if,s)e^P--U{t;s)[PP-^,cv%J 


Since 


we conclude that 


p,(p)-f^{p) = -tvPP-^{u;^,s-io%,s) 


N 


(4.14) 


tr e^P'^{ujN,t-(x%t) 


< 


tiU* {t;0)e^^'^U{t;0){u}N — 0J%) 


1 


+ ^ / ds I dp\V{p)\\trPP--iujN,s-uj%^s)\ tr[W{f,s)e^P-^U{t-,s),e^nco%, 


Ns _ 

and therefore, using the assumption (|2.13p . that 


®^P 1 , I I 
pgR3 1 + \p\ 


tvPP--{uN,t-u:%^t) 
1 


< sup 


pGR3 1 + \P\ 

C 


+ 


/ 


tiU* {t;0)e'^^'^U{t;0){uj]s[ — oj%) 

1 


ds sup 


tr 


Nej, ^~,-3(1 + N)(1 + |p|) 


W{t-,s)PP-^U{t-s),P'P-^ uj%, 


(4.15) 


X sup 


1 + 


To bound the second term on the r.h.s. of (|4.15l) we shall use the following lemma, whose proof is 
deferred to Sect. [H 

Lemma 4.2. Assume that \2.18]} holds true. LetU{t;s) be the unitary evolution generated by the 
Hartree Hamiltonian h{t) = —+ (1/ * p^). There exists a constant C > 0 such that 

sup^ |tr [e*^'^,ZY*(t;s)e“-P+^^-'?Z^(t;s)] u;| < e{\p\ + |g|)e^l*-"l 

uj,r \'^\ 
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for all p,q G Here, the supremum is taken over r G and over all trace class operators uj on 
L^(M^) withtT\u}\ < 1. 

It follows from Lemma 14.21 and from tr|a;^^| = N that: 


tr 


U*{t-s)e^P-^U{t-s),e^P-^ 



< CN£\p\\^ 


(4.16) 


To bound the first term on the r.h.s. of (jd.lSp . we proceed as follows. We choose a function x< € 
with x<{x) = 1 for Ixj < 1 and x<{x) = 0 for |x| > 2. We set x> = 1 ~X<- For an arbitrary 
i? > 1, we decompose 

iTU*{t-,t))H'P'^U{t-,Q){ojN -oj%) 

= ivW{t-, t))HP-^U{t-, 0)x<(-e^A/ii)(a;7v - 

+ tri^*(t;0)e*P'^^^(t;0)x>(-e^A/i?)(wAr - a;^)x<(-e^A/i?) (4.17) 

+ tri7*(t;0)e*P-"^7(t;0)x>(-e2A/i?)(u;,v - a;^)x>(-e'A/i7) 

= I + II + III 


To estimate the last term, we observe that 

-2 A . I / Ji \\, k 


|III| < trx>(-e^A/ii)a;Ar + trx>(-e^A/i?)w^ 

tr (—£:^A)ca^ T tr (—£:^A)ca^ 


< 


R 

N 

'r 


J dxdvv'^Wisi{x,v) + J dxdvv‘^W^{x,v) 


< 


CN 


(4.18) 


from the assumption sup^y ||lLAr||jy 2 < oo, and using that x>(“£^A/i2) < (—e^A/i?). Next, let us 
consider the first term on the r.h.s. of (j4.17p . We write 

I = tr^7*(t; R)e''^'^U{t] 0)x<(— e^A/i?)(l + x^)“^(l + x‘^){ojn — oj%) 


and we decompose 


[(1 + x^){u}N - UJ%)] {x; y) = iV(l + x^) J 


dv 


Wn 


X + y 


v] -W] 


N 


x + y 


(re —Tj) 


= Di{x] y) + D 2 {x-, y) + D^ix] y) 


where 


Di{x-,y) = N 


1 + 


x + y 




ix — v') 


is the Weyl quantization of the function (1 + x^){Wn{x,v) — W^{x,v)) defined on phase-space, while 


D2{x;y) = 


Ne^ (X — y 

~T~ 

Ne^ 


I 


dv 


Wn 


X + y 


,v]-W, 


kfx + y , 

2 ’ 


N 


(x — v) 

e™- . 


/ 


dv 


A,Wn 


x + y 


V I ~ ^vWn 


k(x + y 


(x — v) 
e™ e 
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is the Weyl quantization of {e^/A){/S.^W n{x,v) — A.yW^{x,v)) and 


Dsix-y) = Ne 
= Ns 


x + y X -y 
x + y 


dv 


Wn 


x + y 


,v -W, 


kfx + y 


N 


,V 


I 


dv 


V^Wjv 


x + y 


,v -V,W, 


kfx + y 


N 


,V 


i^-y) 


(x-y) 


is the Weyl quantization of sx ■ ('VyW]\i(x, v) — VyW^{x, n)). We bound the contributions of the three 
terms Di,D 2 ,D 3 separately. We begin with 


trZ^*(t; 0 )e^P'^U{t-, 0 )x<{-e‘^A/R ){1 + x^y^Di 

= trU*{t-, 0)e^P'^U{t-, 0)x<{-e'^A/R){l + x‘^)-^{l - e2A)-^(l - e'^A)Di 

< 11(1 + x^)-\l - e2A)-i||Hs||(l - £'A)Z1 i||hs 

< CVn \\{1 - e^A)Di\\ns 


where we used that 0 < x<(~£^A/i?) < 1. We have 
[{1 - e‘^A)Di]{x-y) 


= N{l-e^A,) 


1 + 


x + y 


dv 


1^0 


x + y 


,v -W^ 


k(x + y 


,v 


(x — v) 


It is not difficult to see that: 

11(1 - e2A)Zli||Hs < C^/N\\{1 + x^)(l + v^){Wn - W^^)h 

+ C'/NeWWn - W^WhI + C^/Ne‘^\\WN - W^\\h 2 


Therefore 


1 




(4.19) 


iiWit] 0)e^P-^U{t; 0)x<{-e‘^A/R){l + x^)-^Di 
The contribution of D 21 on the other hand, can be controlled by 

trZ^*(t; 0)e^P'^U{t-, 0)x<(-e^A/i2)(l + x‘^)-^D2 

< \\x<i-e^A/R)il + x2)-1||hs||I12||hs < Cs^^VnWWoWh^\\ x<i-e^A/R)il + x2)-1||hs 


where 

\\x<{-e^A/R){l + xYY^s 

= tr (1 + x^)-\l - e2A)-V<(-e^A/ii)(l - e2A)-^(l + 

= tr (1 + x2)-i(l - e2A)-i(l - e^A)\l{-e‘^A/R){l - e^A)-\l + x^^ 

<CR ^\\{1 + x^)-\ 1 -s^A)-Y^s 

< CR^N 
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Hence, we conclude that 


tvU*{t- 0)x<(-e^A/i?)(l + x^)-^D 2 < CNRe^ 

We proceed similarly to bound the contribution of the term D^. We hnd 

trZY*(t;0)e'?’'"^/(t;0)x<(-e'A/i?)(l + x2)-iZ)3| < ||x<(-e2A/i?)(l + x2)-1||hs||I?3||hs 

< CNReWWn - W^WhI 

CNRe 


(4.20) 


< 


y/k 


where in the last step we used (I4.2D . The last equation, combined with (14.191) . (|4.20p implies that 


HI < CN 


-|- 6 “h Rs 


Re 


.Vk Vk / 

Analogously, one can show that the same estimate holds for the term II on the r.h.s. of (14.171) as well 
(in this case, we introduce the identity (1 + x^)(l + x^)~^ on the right of the difference ujn — and 
we use the cyclicity of the trace). With (j4.18p , we conclude that 


0)e*P'^Z^(t; 0)(t<;w - < CN 

Choosing R = s~^, we obtain 

tTW{t-,0)e^P-^U{t;0){ujN - uj%) 


1 „ 2 1 
— 1 = + £ + Re H— -j=z + — 

Vk Vk R 




Inserting this bound and (14.161) in (j4.15l) and applying Gronwall’s lemma, we obtain 

1 


sup 
pgk3 1 + IpI 

which concludes the proof of (I4.13|) . 




(4.21) 

□ 


Comparison of with ww,*- We now compare the Vlasov evolution of the regularized initial 
data with the Vlasov evolution of the original data. We claim that there exists a constant C > 0, 
depending on sup^ ||ITAr ||//2 but not on the higher Sobolev norms, such that: 


CeCA 


IIAat,* - A^,il|HS = - lV^,t||2 < 


To prove this, let 


/ 


Pt{x) = / dvWN,tV,v) and pW = / dvWVi-{x,v) 


I' 


(4.22) 


(4.23) 


be the densities associated with WN,t and For t e M, we denote by {Xt{x,v),Vt{x,v)) and by 

{X^{x,v),V^{x,v)) the flows satisfying the differential equations 


Xt(x,v) =2Vt{x,v) 

Vtix,v) =-V{V*pt){Xtix,v)) 


(4.24) 
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and 


r XHx,v) =2 VHx,v) 

\v,Hx,v) =-V{V*p^)iXHx,v)) ^ 

with initial data given by, respectively, Xo{x, v) = Xq{x, v) = x, Vo(x, v) = Vq{x, v) = v. We compare 
the two flows (W) Vt) and {X^^ Vt)- We have 

j^{Xt- X^){x,v) = 2{Vt-Vl^){x,v) 

j^{Vt - V,^)ix, v) = -X{V * pt){Xt{x, V)) + V(y * v)) 

and therefore 


l(Xt-X,^)(x,v) 

±{Vt-V>^)ix,v) 


< 2 


Vt{x,v) - Vt^{x,v) 


<C\\pt-p^\\i + C Xt(x,v)-X,^{x,v) 


where we used the assumption (12.131) . Gronwall’s lemma implies that 


Xt(x,v) - X^(x,v) + Vt{x,v) -V^{x,v) 


< Ce^\^\ 


[ ds\\p. 
Jo 


S Psill 


<( 76 ^ 1*1 [ ds\\Ws 

Jo 


(4.26) 


- w: 


We will also need to control the difference between derivatives of the flows {Xt{x,v),Vt{x,v)) and 
{X^{x,v),V^{x,v)). Integrating the flow equations (j4.24l) . (I4.25p . we have 


XxXt{x,v) = 1 + 2 


/ XxVs{x,v)ds 

Jo 


V, 


,Vt{x,v) = - f X^{V *ps){Xs{x,v)) ■X,Xs{x,v)ds 
Jo 


(4.27) 


which implies that 

\S/xXt{x,v)\<l + 2 [ ds\VxVsix,v)\ 
Jo 

\^xVt{x,v)\<C f ds\V^Xs{x,v)\ 

Jo 

and hence, by Gronwall’s lemma, that 

\V^Xtix,v)\ + \V,Vt{x,v)\ <e^l‘l 

Analogously, we also find 

|V,W(x,u)| + |V,17(x,u)l <e^l*l 

and 


\X^X^{x,v)\ + \X^VtHx,v)\<e^\^\ 
|V,A,"(x,u)| + |V,G,'=(x,u)| <e^l*l 


(4.28) 

(4.29) 
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Moreover, from (I4.27p . we obtain 


V,Xt(x,v)-V,Xt^(x,v) < 2 fds XMx,v) -^xV,Hx,v) 

Jo 


and thus 

V,Vt{x,v)-V,V’^{x,v) 


ft 


< ds 

Jo 

rt 

< _ 

'0 


X^{V *ps){Xs{x,v)) ■X,Xsix,v) - X^{V *r,)iX^{x,v)) ■X,X^ix,v) 

c[ ds\\ps-'^\\i + C f (is|x^(a;,u)-X^^(x,u)||V3;X^(x,u)| 

Jo Jo 

+ C [ ds\XMx,v)-^xX^{x,v)\ 

Jo 

To get the second inequality, we used that 

||W* p.lloo < ||W||oo||VW^||i < Cb^^^^WWnWhI (4.30) 

Using (I4.28P and (|4.26l) . and applying Gronwall’s lemma, we conclude that 

X^Xtix, v) - (x, u)| + I V,Ut(x, v) - X^V’^ix, v) 

pt pt ps (4.31) 

<Ce^l‘l / ds||p,-p^||i + Ce^l*l / ds dr||p,-p^||i 

Jo Jo Jo 

Similarly, we can also show that 

XMx, v) - XyX^ix, u)| + I V,Ut(x, v) - X,Vt\x, v) 

pt pt ps (4.32) 

< / ds IIps - p^lli + / ds / dr ||p^ - p^||i 

Jo Jo Jo 

Next, we control the norm of the difference Wjsf^t — ^Nf end, we write 

\\WN,t - wyu = I dxdv\WN,t{xx) - W^,tix,v)\ 

= j dxdv\WN{X-t{x,v),V-t{x,v))-W^{X^^{x,v),V^^{x,v))\ 

< J dxdv\WN(,X_t{x,v),V-t{x,v))-W^{X_tix,v),V-t{x,v))\ 

+ j dxdv\W^{X_t(x,v),V-t{x,v))-W^{X’lt{x,v),V\{x,v))\ 

Using that the Vlasov dynamics preserves the volume in phase-space, we get: 

\\WN,t-W^,t\\i < \\WN-W^h 

+ j dxdv ^ dX^W^ (\{X_t{x,v),V.t{x,v)) + {I - \){X\{x,v),V^^{x,v))) 

< lltUv - VU^IIi + J dXdxdv \{XxW^){x{x,v,X),v{x,v,X))\\X_t{x,v) - X^i.{x,v)\ 

+ l(V^fU^) (x(x,u, A),i;(x,u, A)) ||U_t(x,u) - V^t{x,v)\ 
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where we introduced the notation 


x{x, V, A) := XX-t{x, ?;) + (1 — X)X\{x, v) , v{x^ v, A) := XV-t{x, v) + (1 — X)V\{x, v) 
From (14.261) . we obtain 

\\WN,t-W^th < \\Wn-W^\\i + C f dse^^^WN,s-W^s\\i 

Jo 

X J dXdxdv I {XxW%) {x{x, v, X),v{x, v, A)) | + | (V^VF^) {x{x, v, X),v{x, v, A)) 

We observe that 

J dxdv\{VxW^){x{x,v,X),v{x,v,X)) \ = J \{VxW^){x,v)\ dxdv 

with the Jacobian 

'X7 xX-t X xV-t 


(4.33) 


(4.34) 


J = det 


A 


X xX-t XyV-t 


+ (1-A) 


XxX^t XxV\ 
XxX’l, x.vj^t 


To estimate the determinant J{x,v) in (I4.34p . we proceed as follows. For a fixed constant C > 0 
(that later will be chosen large enough), let us define t* > 0 such that: 


^3^ 

Vk 


= 1/2 


We claim that, for all \t\ < t*, 


\\WN,t-Wy\i<C^ 


We prove (14.361) for t > 0 (the case of t < 0 can be handled similarly, of course). We set 


to = inf |i 




and we proceed by contradiction, assuming that to < t*. At time t = 0, we have: 
\\Wn — W^\\i< j dxdvdx'dv'gk{x — x',v — v')\Wn{x,v) — Wn{x',v')\ 

- I — W]sfix, v] 


(4.35) 


(4.36) 


(4.37) 


< 


< 


CO 

J “ 

1 

fd 

CO 

J ^ 

X 

|r 

.Vk 

Stt 

1 




dxdvdrds / dXe 
0 

VxWv ( X + X—j=,v + X—j^ 


+ 


y/k 


VtilF/v ( X + ^ 




(4.38) 
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where in the last line we estimated the L^-norms by proceeding as in (I4.30p . Since, moreover, t —)• 
and t are continuous in the L^-topology, by choosing C = 2C in Eq. (I4.37p . we conclude that 

to > 0. The continuity property is a standard fact (see e.g. HI)- 

By dehnition, for 0 < t < to, we have (I4.36P and therefore, from (I4.3ip and (|4.32l) . 


and 


\V^X_t{x,v)-V,X^^{x,v) \ + \ V^V.t{x,v) -X^V^,{x,v)\ < 
\XyX_t{x,v) - X,X’lt{x,v)\ + \X,V.tix,v) - X,V^tix,v)\ < 


e2C'|4 

'/k 

e2C\t\ 

Vk 


Writing 

J(x, v) = det 
and using that 

we conclude that 


^xX-t XxV-t 
^vX-t XyV-t 


det 


+ (1-A) 


VxXK - VxXt VxV^ - VxV- 


x^-t ~ x^t x^-t ~ v^r-i 


'^xX-t VxV-t 
^vX-t V vV-t 


= 1 


\\J{x,v)\-l\<C^ 


„3C7|t| 

'Jk 


if the constant (7 > 0 is large enough. From (I4.35p . and from the assumption to < t*, we conclude 
that 

|J(x,u)| > 1/2 

for all 0 < t < to- Eq. (14.341) implies: 


I dxdv\{VxW^){x{x,v,X),v{x,v,X))\ <2 J dxdv\VxW^{x,v)\ < C||W^||^i <C\\Wn\\hI 
for all 0 < t < to. Similarly, we obtain 

JdxdvliVvW^) (x(x,v,X),v(x,v,X)) | < C\\Wn\\hI 
Plugging the last two bounds in the r.h.s. of (|4.33l) . we hnd that 

WWN^t-W^NtWi < \\Wn-W%\\i + C f dse^\^\\\WN,s-W%,\\i 

Jo 

for all 0 < t < to- Eq. (I4.38p and Gronwall’s lemma imply that, if the constant (7 > 0 is sufficiently 
large, 

~ e^l*l 

for all 0 < t < to, in contradiction with the definition of tg. This shows that to > t*. Repeating the 
same argument for t < 0, we obtain that 

~ e^l*l 

(4.39) 
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for all |f| < f*. From (I4.26p . we also find that 


Xt{x,v) - X^{x,v)\ + \Vt{x,v) 


V^\x,v)\<C^ 


'/k 


for all \t\ < t*. Moreover, Eqs. (I4.3ip and (|4.32l) imply that 


\J{x,v)\ > 1/2 


(4.40) 


(4.41) 


for all \t\ < t* and for all x,v E M^. 

Finally, we control the difference ™ L^-norm. To this end, we observe that 

WN,t-W^NM = j dxdv\WNiX_t{x,v),V.t{x,v))-W^{X^_,{x,v),V\{x,v))\^ 

<2 j dxdv\WN{X_t{x,v),V.t{x,v)) - W^iX_t{x,v)V-t{x,v))\^ 

+ 2 1 dxdv\wUx_t{x,v),V.t{x,v)) - wUx’lt{x,v),V\{x,v))\^ 
Using that the Vlasov dynamics preserves the phase-space volume, we get, for all |f| < t*: 

\\WN,t-wy\^2<nwN-w^\\l 

-|-2 y dX J dxdv ^\{VxW^) {x{x,v, X),v{x,v, X)) \‘^\X-t{x,v) — X^^{x,v)\‘^ 
+ |(V^1U^) {x{x,v,X),v{x,v,X)) \‘^\V-t{x,v) - V^i-{x,v)\^'^ 


<2||lTw-lU^||i 

C<4g4C|t| rl 


+ 2 

< 10 


k 

(-4g4C|t| 


dX / dxdi}[|(V,lF^)(x,{i)|" + |(V,lT^)(x,i))|' 


\J{x,v)\ 


k 


ll^ivll 


m 


To get the first inequality we used the estimate (I4.40F while to get the last one we used (14.411) . By 
definition of t*, we conclude that, after an appropriate change of the constant C > 0, 


WN,t-W^^\\2 < 


CeClil 

y/k 


for all f € M (recall that the bounds ||lTAry|| 2 , ||lVj^j ||2 < C are trivial, since the Vlasov equation 
preserves the L'p norms). This concludes the proof of (j4.22l) . 

Proof of Theorem 12.21 We have, using (14.5p . ()4.10l) . (14.2211 : 

< CVi/2 exp(C exp(C|f|)) f 1 + ^(e\/fc)^' 


(4.42) 
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for a constant C > 0 that depends on sup^y but not on the higher Sobolev norms. Choosing 

k = we conclude that 

\\^N,t - WAT^tllnS < CN^^'^e exp(C'exp(C'|t|)) 

as claimed. This concludes the proof of Theorem 12.21 □ 

Proof of Theorem 12.31 Let WN,t be the solution of the Vlasov equation with initial data Wn- We 
estimate 

\\WN,t - Wth < \\WN,t - WN,th + \\WN,t - Wth (4.43) 

The first term can be bounded by Theorem 12.21 In particular, (I2.15|) implies that 

\\WN,t - WN,th < C'eexp(Cexp(C|t|)) (4.44) 

As for the second term on the r.h.s. of (j4.43ll . we have to compare two solutions of the Vlasov 
equation, with slightly different initial data. But this is exactly what we did in Step 3 of the proof of 
Theorem 12.21 The only ingredients that we used there were a bound for the and for norm of the 
difference of the initial data. Now, by assumption we have ||lTo — lTv||i < kn,i, \\Wo — WAlb < kn ,2 
and ||Wo||j |^2 < C. Therefore, the arguments used in Step 3 of Section U] imply that 

\\WN,t - Wi||2 < C{kn,i + Kiv,2)e^l*l 
Together with (I4.44|) . we conclude that 


\\WN,t - Wth < C'eexp(C'exp(C'|t|)) + C{kn,i + k 7 V, 2 ) exp((7|t|). 


□ 


5 Convergence for the expectation of semiclassical observables 

Here we prove Theorem 12.41 and Theorem 12.51 To show Theorem 12.41 we make first the additional 
assumption that the Wigner transforms Wn of the fermionic operators un are so that sup^y ||ILAr||^4 < 
oo; later, we will relax this assumption with an approximation argument. 

Case supjv ||WAr||j|^4 < oo. We use the expression (|3.4p for the difference ujN,t — ^N,t to write 


= - f tTe^P-^+'^-^^U{t;s)[V * {ps - p.), Z^*(t; s)ds 

£ JO 

1 /■* 

+ - / iihP-^+P-^^U{t-,s)BsU*{t-s)ds 
£ Jo 


(5.1) 


with Bg as defined in (13.31) . We start by considering the first term on the r.h.s. of (15.11) . We have 
tr U{t- s)[V * {ps - p.), WTv,.] U*{t; s) 

= j dz{ps{z)-ps{z))iTe''P''^^'^''^^U{t]s)[V{x - z), 0 JN,sW*{t]s) 

= j dkV{k) j dze-^’^'hPsiz) 

= ^ J dkV{k)tre~'’^'hcoN,s - u}N,s)tr U{t-, s)[e'^''^,ujN,sW*it; s) 
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Hence 


* (ps - ps),uiN,s]i^it; s)ds 


< _ / dk\Vik)\ tre^’^-^{uN,s-^N,s) 


Ctr \uiN,s\ 1 

^ fceRS (1 + 1^1)^ 


tre*^'^(a;Ar,s - (^N,s) 


(5.3) 


X sup • 
u,k 


1 


tr W(i; s)] u; 


where we used the assumption (|2.18p and where the supremum is taken over all fc £ and all io with 
tr |a;| < 1. From Lemma 021 we obtain 


tre^P-^+i-^^Uit;s)[V * {ps - ps),ujN,s]kl{t-, s)ds 


< 


Ctr \ujN, 


1 




tr e''’^''^{ujN,s - tON,s) 


(5.4) 


Consider now the second term on the r.h.s. of (j5.ip . By the cyclicity of the trace, we find 

tj:e^P-^+'i-^^U{t-,s)BsU*{t-,s) = iTU*{t-,s)e^P-^+‘^-^^U{t-,s)B,. (5.5) 

We recall that the kernel of the operator Bg is 

f x + y' 


Bs{x]y) = 


(H * ps){x) - (H * ps){y) - V(H * Ps) 


\ 2 


(x - y) 




Expanding the parenthesis with the potentials in Fourier integrals, we obtain 

{V*ps){x)-{V*ps){y)-V{V*ps)(^^^^^-{x-y) = J dkU{k) ■ {x - y)'^ 


with U = V * Ps- We write 


gifc.x _ gifc.y ^ ^ ^^gifc.(Ax+(l-A)y )■ {x - y) 

Jo dX Jo 


and hence 


(x-y) = [ dX 


^ik-{Xx+{l-X)y) _ 


ik ■ {x — y) 


fl /■! 


= dX i/2)[fe-(x-y)]2 

Jo Jo 


This implies that 


^ pi pi 

Bs='^ dX{X-l/2) dp dkU{k)kikj 


*,i=i' 


Xi 


Xj, e 


i{liX+{l-ti)/2)k-x~^ ^^i{li{l-X)+{l-ti)/2)k-x 
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Therefore, we can bound the absolute value of the second term on the r.h.s. of (15.11) by 

u dA|A —1/2| J dfj, J dk\U{k)\\k\'^ 


< 


*j=i' 




^ rl pi p 

V / dA|A-l/2| dfi dk\U{k)\\k(^ 

i,j=i Jo Jo J 




< 


CtrlwAT^sl f dA: |17(A:)||A;p sup |tr [xj, s)]]a;| . 

J 


(5.6) 


The supremum on the r.h.s. is taken over all indices i,j G {1,2,3} and all trace class operators ui 
with tr IcjI < 1. This term is controlled thanks to the next lemma, whose proof is deferred to the end 
of the section. 


(5.7) 


Lemma 5.1. Under the same assumptions of Theorem \2.4\ there exists C > 0 such that 
sup |tr [xi [xj , U*{t; s) U{t; s) ]] a;| < Ce^(|p| + 




where the supremum is taken over all i,j G {1,2,3} and over all trace class operators on L^(M^) with 
tr Iwl <1. 


Using \U{k)\ < |U(fc)|, the assumption (|2.18p and (15.71) . we conclude that 

i^UP-^+<i-e^U{t-s)BsW{t-s) < Ctr|w^,,|(|p| + |g|)Ve^l*-"l 

Inserting (15.41) and (15.81) on the r.h.s. of (15.ip . we obtain 
1 


p,qGR3 {\P\ + kl + 1)^ 




tr \ ^N,s\ C|t-s| 


1 


(1 + 


<C [ ds 

~ Jo IV u 

+ C f dsti\ujN,s\£e^J~^^ 

Jo 

<C tds SUP_i_ 

-Jo N 7Ui + \p\ + \Q\r 

+ C f dsti\ujN,s\£e^J~^^ 

Jo 


{uj]y,s - ^N,s) 


(5.8) 


\trUP-+<i--^{u;N,s-J^N,s)\ 


(5.9) 
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Now, we estimate the trace norm of con^s (here we need the additional regularity of the Wigner 
transforms of the initial data assumed at the beginning of the proof). We have 

tr = tr |(1 - + x^)(l - 

< ||(l-e^A)"^(l + x^)“^||Hs||(l + a:^)(l-e^A)wAr,s||HS (5.10) 

< C'/N\\{1 + x^){l - e^A)a;Ar,5||HS 

The operator K = (1 + x^)(l — e^A) wtv.s has the integral kernel 

K{x;y) = N{l + x‘^){l - j dvWN,s{^^—^,v^e~"^'^ 

= N{l + x‘^) J dv Wn,s ( ^ ^ ^ ^ 

+ A(l + x^) J dv v‘^Wn,s( J ^ ^ i 

-e^A^(l + x^) j dv {AyWN,s)(^^^-^,v^e~"'''^ 

+ ieN{l + x^) J dvv ■ VyWN,s(—^^,vje~'''''^ 

Writing 

we conclude that 

4 

\\K\\ns < CVNj2e^\\WN,s\\Hi 

j=0 

The propagation of regularity of the Vlasov equation from Proposition IB.II gives us 

4 

\\K\\ns < CVNc^^^^^e^WWnWhI 

j=0 


and thus, with ()5.inp . 

4 

tr|cjjv,sl < 

Hi 

j=0 

Inserting in (|5.9p and applying Gronwall’s inequality, we hnd 


o,gGM3 (1 + \P\ + \d\y 


tre^P-+i--^{coN,t-d;N,t)\ 

4 4 

Cc[Y,e^\\WN\\Hi\Neexp exp(C|t|)) 


(5.11) 


1=0 


1=0 


This completes the proof of Theorem 12.41 under the additional assumption that ||kVv||_ff4 is 
bounded. 
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Proof of Theorem 12.4L We have to relax the condition sup ||lTAr||j|^4 < oo. To this end, we proceed 
as follows. We set 

W^{x,v) = {Wn * gk){x,v) = j dx'dv'gkix -x',v - v')Wn{x',v') 


with 




and we denote by o;^ the Weyl quantization of W^. We recall from (14.311 . that 


(2vr)= 


uj%{x-,y) = J dwdze 


e VkegVk Ulj^e vfe e vfce 


ix;y) (5.12) 


is a fermionic reduced density with 0 < uj^ < 1 and trw^ = N. In fact, (|5.12ll . together with the 
assumption ()2.19ll . also implies that 


I k \ ^ 

trjcjTV — ^ ^ 

V ^ 


(5.13) 


To see this, we write: 


trlwjv-wil < —^ [ dwdze-^^/^e-^^/^tr 

^ ^ - ( 27 r )3 J 


iw-^ ^-V —f=-V -iw- 


e Vkeg^/k ujj\ie v'fc e ~ 


where 


tr 


iw^ ^-V -^-V -iw-^ 

e Vfceevfc wwc 6 ~ 


< tr| [ujn, e ^ I + tr|e™ v^® — a;Ar| 

<tr|[a;Ar,e '^ ^] | + tr| [ww, e *"'''^®]| 

\^\ I I 


this estimate together with the assumptions (|2.19l) implies that: 


tr cjjv — ^n\ — 


CN 1 


dwdze-^^/^e-^^/‘^[\z\ + \w\] , 


y/k (27r)3 

which proves Eq. ()5.13ll . 

We have for all j = 1,... ,4. Choosing k = (|5.11ll implies that 


D,QeR3 (1 + \p\ + \q\y 


tr — cjAT^i) < CA^eexp(C'exp(C'|t|)) 


On the other hand, proceeding as we did between (|4.14p and (|4.15p (replacing the observable 
with e*p-a:+'?-eV^^ obtain 


sup 


p,, 1 + IpI + kl 


- u;%^t) < Cexp{C\t\)tT\u,N - uj%\ 
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With (15.131) . we conclude that (again with the choice k = e 

1 


sup ■ 


- uj'ht) < CNeexp{C\t\) 


P,^l + \p\+\q\ 

Finally, we observe that 

tr - ^%,t) = N [WnAp, Q) - W^.ip, g)) 

and therefore we estimate 

< CN\\WN,t - W^h 

The L^-stability of the Vlasov equation with respect to perturbation of the initial data has been 
already established in the proof of Theorem 12.21 Following the arguments between (14.231) and (I4.39p 
(using the assumption on the Sobolev norm of the sequence Wn), we obtain 

\\WN,t-W^,th<Ce^^'WN-W^h 

Using again the uniform bound HWatIIwui < C, and the choice k = we find 

\\WN,t -W^,t\\i < 

We conclude that 

1 


sup ■ 


|tre*P'^+'*"^'^(wAr,t - u}N,t)\ < C'Veexp(C'exp(C'|t|)) 


{1 + \p\ + \q\r 

for any sequence of initial densities lon satisfying (I2.19P and whose Wigner transforms Wn satisfy 
WWnWw^,^ < C uniformly in A^. □ 


Proof of Theorem 12.51 We write 


WnAp^q) -^t{p,q) < 


^N^tip^q) - WnAp^q) 


+ 


WnAp^Q) - Wt{p,q) 


where W]\f,t denotes the solution of the Vlasov equation with initial data Wn- From Theorem 12.41 we 
know that 


WnAp^q) - ^nAp^q) 


<Ce{l + \p\ + \q\fe^\^\ 


To conclude the proof of the theorem, we need to compare the solutions WN,t and Wt of the Vlasov 
equation, using the fact that the two initial data are close in L^. As in the approximation argument 
used in the proof of Theorem l2.4[ we make use of the L^-stability of the solution of the Vlasov equation, 
established in Step 3 of the proof of Theorem 12.21 Following the arguments between (|4.23p and (I4.39D , 
we obtain 

\\WN,t-Wt\\i<Ce^^A\WN-Wo\\i 

where the constant C > 0 depends only on ||Wo||wm- This implies that 


\\WN,t - fUilloo < \\WN,t - fUilli < Cknc 


Ulil 


Hence, 


WnAp^ a - Wt{p, q) < C{1 + \p\ + |g|)^(e + KAr)e‘^l*l 
which concludes the proof Theorem 12.51 


□ 
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6 Proof of auxiliary lemmas 

In this section we show Lemma 14.21 and Lemma 15.11 


Proof of Lemma \4.^ We dehne the unitary evolution U{ps) satisfying 

iedt U{t; s) = h{t) U{t- s) 

= {h{t) + 2ie^r ■ V + r'^e^)U{t] s) 

We observe that 

sup|tr =sup tr s)] 0) 

UJ U) 

iiU*{s- 0) s)] U{s- 0)a; 


( 6 . 1 ) 


= sup 

U) 


( 6 . 2 ) 


where the supremum is taken over all trace class operators uj on L^(M^) with tr |a;| < 1 and where we 
used the fact that tr \U{s-, 0)wiY*(s; 0)| < tr |ti;|. For a fixed uj and for fixed f G M, we compute now the 
time-derivative of 

iedsiTU*{s- 0) s)] U{s- 0)a; 

= - ixW{s] 0) [h{s), s) s)]] Uis] 0)a; 

- 2e2 tiU*{s-, 0) ir • V s) s)] U{s-, 0)w 

- r2 trf^*(s; 0) ,U*{t; s) s)] U{s] 0)a; 

+ irW{s] 0) [e*"■^ [h{s)M*{t] s) s)]]Z^(s; 0)a; 


Using the properties of commutators, we find 

iedsiiWis; 0) [e*’'•^^Y*(^; s)] U{s] 0)a; 

= - 2e2 trW*(s; 0) ir • V s) s)] U{s; 0)w 

- r2 trf^*(s; 0) s) s)] U{s] 0)a; 

+ trW{s] 0) [W{p s) s), [/i(s), U{s-, 0)a; 


We have 

[/i(s), = {-2ie^r ■ V - . 

Inserting this expression in (16.311 . we get 

iedstiWis; 0) [e*’■•^iY*(^; s)P^-P+^^ PU{t; s)] U{s-, 0)a; 

= 2etiU*{s-, 0) [U*{t; s) s),ir ■ eV] e^''-^U{s- 0)w 

Integrating this equation from time s to time t, we find 


trZ^*(s; 0) s)] U{s; 0)a; 

= trf^*(f; 0) U{t- 0)a; 

+ 2i^ drtrW*(r;0) [Z^*(f; r) r), fr • eV] e*’''^Z^(r; 0)t<; 


(6.4) 
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which implies that 


iiW{s- 0) s)] U{s- 0)w 




+ 2 y dr 


tr \U*{t- T)e^^-P+^^-m{t-T),ir • eV] e^^'^Uir-, 0)ojU*{t; 0) 


Since 


g*a:-p+eV-q] _ ^^-ier-ql2 _ ^ier-q/2'j^ix-(p+r)+e'V-q 


we conclude that, for any trace class operator uj on with tr |a;| < 1, we have 


ft 

- 

/ dr sup 
Is 'xi 

tr 


W{t; T)e^^-P+^^-'^U{t-, T),i^ ■ sV 

r 


UJ 


where, on the r.h.s., the supremum is taken over all trace class uj with tr |a;| < 1. From ([621) , we 
obtain 


sup 

u),r 


tr [e^^-^,W{t-,s)e^^-P+^^-m{t-s)] uj 


<e\q\ + 2 f 

J S 


dr sup 


tr 


\r\ 


UJ 


(6.5) 


Next, we bound the supremum on the r.h.s. of the last equation. To this end, we observe that 


sup 


tr 


■ eV,U*{t-, s) 


UJ 


= sup 


= sup 

U) 


tr 


U*{t- s), • eV 


trZ^*(s; 0) 


W{t- s),i— • eV 

r 


U{s]Q)ujW{s-,Q) 
U{s; 0)a; 


( 6 . 6 ) 


We compute 

i£dstTU*{s; 0) s),ir ■ eV] U{s-, 0)a; 

= — trZ^*(s; 0) [/i(s), [U*{t] s) U{t] s), ier ■'V]]U{s;0)uj 

+ trZ^*(s; 0) [[/i(s), U*{t] s) s)] , ier ■ 'S/]U{s-, 0)a; 

The Jacobi identity implies that 

iedsi^U*{s] 0) s), ir ■ eV] U{s] 0)t<; 

= - iiW{s- 0) pl*{t- s) s ), [h{s ), ier • V]] U{s-, 0)a;. 

We have 


(6.7) 


[h{s),ir ■ eV] = ier ■V{V * Ps){x) = ier 


j dkkV{k)%{k)e^^-^ 
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Hence 


iedgi^U*[s] 0) s), ir • eV] U{s] 0)w 

= -ie- j dkr- kV{k)pt{k)iYW{s]^) ^^(s; 0)a;. 

Integrating from time s to time t, we find 

trZ//*(s; 0) s),ir ■ eV] U{s; 0)a; 

= trZ^*(t; 0) ir ■ eV] U{t-, 0)a; 

+ J dkr ■ kV{k)Mk)tTU*{T-0)[U*{t;T)e^^-P+^^-m{t-T), e^'^-^]U{T-0)uj 

Since 


ix-p+eS/-q 


[e“- 


,ir ■ e 


V] = £r ■ pe“'- 


ix-p+eV-q 


we conclude that, for any trace class operator u with tr |u;| < 1, 


trZ^*(s;0) 


s)e^^-P+^^-'^U{t- s), ■ eV 

r 


1 


U{s] 0)ti; 


= e|p| + / dr sup— tr^^*(T; 0) r) r), ZY(r; 0)a; / dk\k\‘^\V{k)\ 

Js uj,k \^\ J 


From (|6.6|) . we find 
tr 


sup 

uj.r 


U 


i— ■ eV,U*{t-, s) 

\r\ 

ft 1 

<e|p|+C7 / dr sup— trZY*(T; 0) r) r), ^^(r; 0)a; 


Combining this bound with (16.5|) and applying Gronwall, we obtain 


sup 

a;,r 


tr [e^^-^,U*{t]s)e^^-P+^^-'^U{t-,s)] co 


sup 

LO,r 


tr 


• eV,U*{t; s)e^^-P+^^-PU{t; s) 


OJ 


<eN{\p\ + \q\)e^\^-^\ 


□ 


Proof of Lemma 15.fl We observe, first of all, that 


sup 


tr [xi [xj , U* (t; s) U{t]s)]]uj 

= sup |tr [xi [xj , U*{t; s) U{t-, s) f\U{s-, Q)ojU*{s\ 0)| 

UJ 

= sup |tr^^*(s; 0) [xi [xj , U*{t, s) U ft', s) ]]U{s] 0)a;| 


( 6 . 8 ) 
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We consider now the derivative 


iedstrU* [s] 0) [xi [xj , U*{t] s) U{t] s) ]]U{s] 0)a; 

= — iiU*{s; 0) [/i(s), [xi [xj , U*{t; s) U{t; s) ]]] Z^(s; 0)a; 

+ trZ^*(s; 0) [xi [xj , [h{s),U*{t] s) U{t] s) ]]] ^(s; 0 )li; 

The Jacobi identity implies that 

i£dstTU*{s-, 0) [xi [xj , s) U{t; s) \\U{s; 0)a; 

= trZJ*(s; 0) [xi [[xj, h{s)],U*{t] s) U{t] s) ]] ZJ(s; 0)a; 

+ trZJ*(s; 0) [[xj, h{s)], [xj,U*{t; s) U{t, s) ]]U{s-, 0)a; 


Since [xj,h{s)] 


(and since [V^^.,Xj] = Sij is a number), we conclude that 

iedstTU*{s-, 0) [xj, [xj , U*{t; s) ZJ(t; s) \\hl{s; 0)a; 

= etrZJ*(s; 0) [eV^^., [xi,U*{t; s) U{t; s) \\U{s; 0)a; 

+ etiU*{s] 0) [eVxi, [xj,U*{t-, s) U{t] s) ]] ^(s; 0)ci; 


Integrating over time, we find 

tr^J*(s; 0)[xj, [xj,U*{t-, s)]]^J(s; 0)a; 

= trZJ*(t;0)[xi, [xjyP-^+'i-^^]]U{t-Q)uj 

+ iy drtrZJ*(r; 0)[eV2;^., [xj,ZJ*(t; r)]]^J(r; 0 )a; 

+ i f dTtTU*{T;0)[£Vxi,[xj,U*{t-,T)e'^^'^~^‘^'^^U{t;T)]]U{T-,0)u} 

J S 


Since 

we find 


[x,, [x,-,e'^’-"+'^-"^]] = 


|trZJ*(s; 0)[xj, [xj,U*{t; s)]]ZJ(s; 0)a;| 


< e 


'|gp+ [ dr sup |tr [eVo;^., [xj,ZJ*(t;r)e*^'^+'^''^'^ZJ(t;r)]]u;| 

Js LU.iJ 


for all trace class to with tr |a;| < 1. From (|6.8p . we obtain 
sup |tr [xi, [xj^U*{t] 5)]]a;| 

f dr sup|tr[eV 3 ;^.,[xi,iY*(t;r)e*^’'^+'''^'^ZJ(t;T)]]a; 

Js 

where the suprema are taken over all trace class uj on with tr |a;| < 1. 

Next, we look for an estimate for 

sup |tr [eVj;^., [xi,X(t;s)e*^'^+'?''^'^ZJ(t;s)]]w| 

= sup |trZJ*(s;0)[£Va;., [xj,ZJ*(f; s)]]ZJ(s; 0)ti;| 


(6.9) 
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To this end, we compute the derivative 

ie5strZ^*(s;0) [eV^;;, [xj,U*{t-,s)e''^'^^^'^'^U{t-,s)]\U{s-,Q)oj 

= — trZY*(s; 0) [/i(s), [eV^i, [xj,U*[t] s) U{t] s)]]]U{s]Q)ijj 

+ trZ^*(s; 0) [eVa;., [xj, [h{s),U*{t] s) U{t] s) ]]]Z^(s; 0)a; 

= tr^^*(s; 0) [[eVa;., h{s)], [xj,U*{t; s) U{t] s) ]] ^(s; 0)u; 

+ tr U* (s; 0) [eVa:^, [[xj, /i(s)], U* (t; s) ZY(t; s) ]] Z^(s; 0)a; 


With 


we obtain 


[eVa;;,/i(s)] =eVxiiV * Ps)ix) = e J dkkiV{k)psik)e'‘'^''" 


( 6 . 10 ) 


i£dsticU*{s; 0) [e'Vxi, [xj,U*{t; s) U{t; s) ]]Z^(s; 0)a; 

: J dkkiV{k)ps{k)ivW{s]^) [xjM*{t-,s)(^P'^+‘i-^^U{t-s)]]U{s;Q) 


= e 


+ eiTU*{s-, 0) [eVa:,, [eVx,M*{t] s) s) ]]U{s- 0)a; 


Using the identity 

= C dX-^ ^i\k-x j^^i{i-\)k-x ^ dXe^^^-^[ik ■ x, (6.11) 

Jo dX Jq 

we conclude that 

i£dstTU*{s-, 0) [eVxi, [xj,U*{t; s) U{t] s) ]] ^(s; 0)a; 

= e J\x j dkkikiV{k)psik)tTU*is-,0)e^^’^-^[xi, [xj,U*{t]s) U{t-,s)]]e^^^-^^’^-^U{s;0)io 

+ etrZ^*(s; 0) [eVa,,, [eV^„iY*(t; s) U{t; s) ]]Uis; 0)w 

and hence, after integrating over time. 


tr^^*(s;0) [eVxi, [xj^U*{t-,s)e'’^'^^‘^'^"^U{t-,s)]\U{s-,^)ijj 


< 


[eVx,, [xj, 


+ 


+ 


[ dr [ dX [ dk\k\‘^\V{k)\ 

Js Jo J 

X trZ^*(T; 0) [xj,U*{t-, r) r) ]] 0)w 

f dr iiU*{T-,0)[£Vxi, [£Vxj,U*{t;T)e^P''"~^'^'^^U{t]T)]]U{T]0)Lo 

J S 


Since 


[eVa.,, [x„e*P'"+'?-^]] = 
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this implies that 


sup 


tr [eV,,, s) s) ]]a; 


< e"^|p||g| + C / dr sup \tr [xi, [xj,U*{t-,T) U{t-,T)]]oj\ 

Js 


( 6 . 12 ) 


+ dr sup 

JS Uj,i,j 

Finally, we need an estimate for 


tr [sVxi, [eV^. r) U{t;T)]]oj 


sup 


tr [eV,,,, [EV,.,Wif, s) s) ]] cu 


= sup 

^,hj 


tiW{s, 0)[eV,,, [eV,,,U*it- s) U{t- s) ]]U{s- 0) a; 


(6.13) 


Hence, we compute the derivative 

iedsixU*{s- 0)[eV^,, [eV^.,Z^*(t; s) s) ]]U{s- 0) w 

= - tiW{s- 0)[h(s), [eV,,, [eV,,,^/*(t; s) s) ]]]U{s- 0) w 

+ trZ^*(s; 0)[eV,,, [eV,., [h{s)M*{t- s) s) ]]]Z^(s; 0) w 

= tr W {s- 0) [[eV,,, h(s)], [eV,,, W (t; s) Z^(t; s)]]^/(s; 0)w 

+ trZ^*(s; 0)[eV,,, [[eV^.,h{s)],W{t- s) U{t; s)]]U{s-, 0)a; 

From (jh.lOll . we find 

iedstr U* (s; 0) [eV^^, [eV^j , U* {t; s) U(t; s) ]]U{s]0)uj 

= e y dkkiV{k)ps{k) iTU*{s] 0)[e*^'^, [eVx^M*{k, s) U{t, s)]]ZY(s; 0)w 

+ e y dkkjV{k)ps{k) tiU*{s-, 0 )[eV 3 ;., [e*^'^,Z^*(t; s) U{t] s)]]iY(s; 0)ti; 

In the first term on the r.h.s. of the last equation we use (|6.1ip . In the second term, on the other 
hand, we notice that 

tTU*{s; 0)[eV,,,,[e*^■^^^*(^; s) U{t- s)]]ZY(s; 0)w 

= trZ^*(s; 0)[e*^■^ [eVx,M*{t] s) s)]]Z^(s; 0)a; 

+ iiU* (s; 0) [[eV,,,, , U* {t; s) U{t-, s)]U{s; 0)w 

Again, the first term on the r.h.s. of the last equation can be handled with (|6.1ip . As for the second 
term, we use that [eVa,., = iekie'’^'^. Integrating (I6.13p over time, we find 

tiW{s; 0)[eV,„ [eVx,M*{t- s) U{t, s) ]]U{s; 0) w 

< |tr^^*(t;0)[eVa.„ [eVa.,, ]]Z^(t;0)a;| 

+ C f dr sup |tr [eVa;j, [a:j,^^*(t; r) r)]]a;| 


's 


+ Ce 


dk\V{k)\\Mk)\\kfdk 


1 


dr sup^ tr[e*^-^ZY*(t;r)e*P'^+'^-"^ZY(t;r)]w 
Is oj,k \k\ 
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To bound the integral involving the potential in the last term on the r.h.s. of the last equation, we 
use (|2.18p with ||yOs||cx) < 1- From 

and from Lemma iia we obtain 


iTU*{s- 0)[eV,,, s) s) ]]U{s-, 0) a; 


<C£‘^ + C f drsup|tr [£Vxi,[xj,U*{t;T) U{t-,T)]]u}\ 

Js uJ,iJ 


(6.14) 


sup 


Combining (16.91) and (16.121) with the last equation and applying Gronwall lemma, we deduce that 

tr [xi, [xj,Wit-s)e^P-^+^-^^Uit;s)]]uj\ < C£^i\p\ + 
tr [eV,,, [xj,U*{t;s)e^P-^+'^-^^U{t;s)]]u;\<Ce\\p\ + \q\fe^^^^^^ 
tr [eV,„ [eVx„U*{t-, s) U {t- s) ]] o; < Ce\\p\ + 


sup 


sup 


□ 


A Well-posedness of the Vlasov equation for signed measures 

The goal of this appendix is to show that the arguments of [8] can be extended to prove the well- 
posedness of the Vlasov equation (11.101) for initial data given by signed measures. 

Following the notation of [8], let At denote the space of all hnite signed measures on the Borel 
(7-algebra For an open interval A C M, we denote by A4 a the set of all families M = {pt}teA, 

with /if £ At for all t £ A such that, for all bounded intervals A' C A, there exists Ca' > 0 with 
supfgA' ll/^tll < C'a', and such that the function 


(VV * /if) (x) 


W(x — x') dpt{x', v') 


is continuous in t £ A, for all x £ (V denotes the interaction potential entering the Vlasov 
equation ()1. 101) 1. For C > 0, we also denote by the set of families M = {pt}teA with 

||/if|| = supsgg(K6) \p{B)\ = K for all t £ A. 

Defining A, i? : x ^ x by A{x,v) = {2v,0) and B{x,v) = (0, W(x)) and, for every 

/i £ At, 


Bfj,{x,v) = / B{x — x',v — v')dp{x',v') 


we say that a family M = {pt}teA £ A^a is a weak solution of the Vlasov equation on the interval A 
if, for every test function h £ (D(R®) in the Schwarz space. 



/i(x, v)dpt{x, v) 
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is differentiable in t E A and 

^/Xf(/i) = utiiA + B^^)Vh) 

It is easy to check that, if the weak solution /r* has a density Wt{x,v), differentiable in t, then Wt is 
a solution of the standard Vlasov equation (jl.lOp . 

Proposition A.l. Let V € For any finite signed measure G Ai, and every open interval 

A C M with 0 G A, there exists a unique weak solution M = {fat}t£A of the Vlasov equation on A 
with p,t=o = ■ 

Proof. We follow the strategy of [8], adapting it to the case of signed pfi. We will use the variable 
z = {x,v) G M®. For M = {pt}t^A £ -A^A) we define 

z) = A{z) + B^fiz) (A.l) 

and we consider the solution of Newton’s equation 

^z{t) = GM{t,z{t)) (A.2) 

We denote by ZM{t,u) the solution of ()A.2h . with initial data zm{^-,u) = u. 

For a fixed G Ad, we dehne the map T : A^a —^ A^a by 

{TM)t{E) = / ({n G : ZM{t,u) G E}) 

for all E G As in [H], it is easy to check that M G Ad a is a weak solution of the Vlasov equation 

with initial data p^ if and only if M is a fixed point of T, i.e. if TM = M. 

Hence, to prove Proposition lA.ll it is enough to show that T is a contraction on AdA- In fact, 

since clearly TM G AdA(||Ai'^||), for all M G Ad a, it is enough to show that the restriction of T to 

AdAdll^^ll) is a contraction, with respect to an appropriate metric, that we are now going to define. 
For two signed measures /i,/z' G Ad, we define 

d{ti,l^') = dKR{l^+,l^'+) + dKR{p-,P-) (A.3) 

where p = p^ — p- is the Jordan decomposition of p in its positive and negative parts and where dxR 
is the Kantorovich-Rubinshtein metric, defined by 

dKR{R,R')= inf p{zi, Z 2 )dm{zi, Z 2 ) 

7nGN{i',u') J 

where p{zi,Z 2 ) = min(| 2 i — Z2\,l) nnd N(n,n') is the space of all positive measures m on 
such that m{E x M®) = ^{E) and m(M® x E) = v'{E) for all E G BfMfi). Furthermore, for M = 
{pt}t€A,M' = {p[}teA G Ma, we define 

d{M,M')= [ d{pt,p[) 

Ja 

It is easy to check that ()A.3p defines a metric on Ad (A). 

We claim that, for |A| small enough, 

d{TM,TM') <^d{M,M') (A.4) 
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for all M,M' G A^A(llAi'^ll)- To prove (IA.4P we observe that, for all G M., 

B^{z) — B^i{z) = J B{z — w)d^{w) — J B{z — 'w)d^'{w) 

< J {B{z - wi) - B{z - W2))dm+{wi,'W2) 

— J {B{z — wi) — B{z — W 2 ))dm-{wi,W 2 ) 

for any m+ G G A^(/r_. Recalling that B{x,v) = (0, VR(x)) and the assumption 


V G we find 

\B^{z) - B^>{z)\ < j\B{z - wi) - B{z - W 2 )\dm+{wi,W 2 ) 

+ J \B{z — wi) — B{z — W 2 )\dm-{wi,W 2 ) 

< Cy p{wi,W2)dm+{wi,W2) + J p{wi,W2)dm-{wi,W2) 

Since the inequality holds for every m+ G N{p^,p'^) and m_ G N{p-,p'_), we conclude that 

\Bi,{z) - Bf,>{z)\ < Cd{p,p') (A.5) 

for all p, p' & M and all z G M®. 

Furthermore, recalling the definition (jA.ip . we observe that there is a constant C, depending on 
||/i®||, such that 

\Gm{z)-Gm{z')\<C\z-z'\ (A.6) 

for all z,z' G M® and for all M G A4a(||m°||)- 

For M,M' G A4a(||/^°||) and n G M®, we dehne the quantity 


a{M, M', u) = sup |zm(T y) — ZM'{t, ^i)| 
t€A 


With ()A.6p . we obtain 


\zM{t,u) - ZM'{t,u)\ < / \Gm{s,zm{s,u) - GM'{s,ZM'{s,u))\ds 

Jo 

< / \GM{s,ZMis,u) - GM{s,ZM'{s,u))\ds 

Ja 

+ / Gm{s,zm'{s,u) - GM'is,ZM'{s,u))\ds 

Ja 

for all t G A. Combining (jA.Sp and (lA.Gh . we find 

\zM{t,u) - ZM'{t,u)\ < G / \zm{s,u) - ZM'{s,u)\ds + Cd{M,M') 

Ja 

Taking the supremum over t, we conclude that, for sufficiently small |A|, 


a{M,M',u) < 


G 


1-C|A 


■d{M, M') 
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We are now ready to bound d{TM,TM'). For M,M' G A^A(||Ai°||), we notice that 

{TMt)±{E) = 4 ({n G M® : ZM{t, u) G E}) 
for every E G Now, let 

mt,±{E) = [{u G M® : {zM{t,u), ZM'{t,u)) G F}) 

for every F G ;S(]R^^). Then we have mt^± G N{{TMt)±, {TM[)±) and 


p{zi,Z2)dmt,±{zi,Z2) = J p{zM{t,u),ZM'{t,u))dp^{u) 

< J a{M, M',u)dp^{u) 

rd{M, M') 


c\\pI\ 

1-C|A|' 


This implies that 

dKR{{TM,)^,{TMi)^) < 

and therefore that ^ 

d{TM,TM') < 

1 - (7|A| 

Hence, for |A| sufficiently small, we obtain (IA.4p for all M,M' G A4a(||/^®||)- This proves that T 
defines a contraction on A4 a(||/^®||) and implies the existence and the uniqueness of a weak solution of 
the Vlasov equation, for |A| sufficiently small. The argument can then be iterated to obtain existence 
and uniqueness for all times. □ 


B Regularity estimates for solutions of the Vlasov equation 

In the next proposition we estimate the weighted Sobolev norms of the solution at time t of 

the Vlasov equation in terms of their value at t = 0. 

Proposition B.l. Assume that 

j dp\V{p)\{l + \p\‘^) <oo (B.l) 

Let Wt be the solution of the Vlasov equation M.KA) with initial data Wq. For k = 1,2,3,4,5, there 
exists a constant C > 0, which depends on ||lTo||j:^| but not on the higher Sobolev norms, such that 

\\Wt\\Hi<Ce^''Wo\\Hl (B.2) 

Proof. We use a standard argument. We denote by <I>t(x,u) := {Xt{x,v),Vt{x,v)) the solution of 
Newton’s equations 

Xt{x,v) = 2Vt{x,v) 

Vt{x,v) = -X {V *pt) (Xt{x,v)) 
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with initial data Xq{x,v) = x and Vo{x,v) = v. Here pt{x) = f dv Wt(x,v). We can rewrite Newton’s 
equation in integral form 


Xt(x,v)=x + 2 / dsVs(x,v) 

Jo 

Vt{x,v)=v-[ V{V * ps){Xs{x,v)) 
Jo 


(B.3) 


In the following, it will be convenient to introduce the following shorthand notation; 

ll^i^^lloo := max||V“W(x,u)|| 
l«l=j 

||V'/'^^||oo := max ||V“Vt(x,u 
l«l=j 

U) II I IIT/(i) I 


(B.4) 


||#||oo := ||x(^'^||oo + ||Vi 


t Iloo • 


( 7)11 

In general, to control it is sufficient to control ||<I)j ^|oo for j < k. In fact, it is not difficult 

to see that: 


\\Wi 


tWnk 


f dxdv {I + x^ +u^)^|V"Wo(X_t(x,7;), V-t{x,v))^ 

\a.\=k'' 

E / dxdv {I + x^ X v^)‘^\{X^WQ){X_t{x,v), V-t{x,v))\ 

1 / 01^7 rxi rv I /I I ^ 


\P\<k «l,-",a|/J| 

ai+...+a^p\=k, |oi|>l 


< 


X I (X_i(x, v), V.t{x, v)) I' • • • I V“l/5l {X_t{x, v), V.t{x, v)) f 


(B.5) 


ri=l 

miH-...H-mn=Ai, mi>l 


to get the last step we performed a change of variables and we used that, by Gronwall’s lemma together 
with (IB..Sp and ||VH||oo < 00 : 

1 + X^{x, v) + V^{x, v) < + v^) (B.6) 

We start by estimating ||IW||j|^i. To this end, we need to control For any multi-index a with 

|q:| = 1, we obtain from ()B.3p that 

||V“W||oo <1 + 2 f ds||V“H,||oo 
Jo 

||V“Ht||oo < 1+ [ ds\\X\V*ps)oXs-X^X,\\^ 

Jo 

<1 + C [ ds||V“X,||oo 

Jo 

where we used that ||V^(H * Ps)||oo < ||V^H||oo||/Os||i, and ||ps||i < ||Ho||i < C'HHoHji^o (see (|4.3Up b 
Gronwall’s lemma, together with the assumption ||V^H||oo < 00 , implies that 

ll^r^lloo < (B.7) 
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where the constant C depends on ||Wo||j|^o, but not on the higher Sobolev norms. Thanks to ()B.5I) . 
the bound ()B.7p immediately implies 


Wt\\l^<Ce^\^\\\Wo\\li 


(B.8) 


where the constant C depends on ||kTo||j:/o, but not on the higher Sobolev norms. This concludes 
the proof of (IB.2h with k = 1. Next, let k = 2. As before, we start by considering the derivatives 
V"Xt, V“Vt, now for |a| = 2. We have: 


|v“(v(y * ps) o x,)||oo < I|v^(y * Ps)\U\xi^^ IlL + * P.)llool|x?^lie 

<C'e^W||Wo||m + l|V2(y*p,)|U||X(2)||^ 


(B.9) 


where in the last step we used ||V^(y */9<j)||oo < ||V^17||oo||V/?s||i < ||lTo||j:^i, and we estimated 

||V/Os||i < C'll^sll/fi < ||l^o||_ffi• This, together with the estimate ()B.7D . implies: 


l^f^lloo < 2 


I Iloo ^ 


[ rf^lK^'^lloo 
Jo 

C fds\\XP\\^ + Ce^\^\\\Wo\\Hl 
Jo 


thus, by Gronwall’s lemma: 


||<f>f^||oo<t^e^'‘'||W^o|lHi 

Therefore, proceeding as in (IB.Sp . we get 

\\Wt\\H^<Ce^\^\\\Wo\\m 


(B.IO) 


(B.ll) 


where the constant (7 > 0 is allowed to depend on HWoH/^i, but not on the higher Sobolev norms. 
This concludes the proof of ()B.2p for k = 2. Consider now k = 3,4,5. We will use that, for |a| = k: 


|V“V(y*p,)||oo<C||VV||oo liv^p, 


'sill 


1=fc-l 


(B.12) 


^CllvVlIoolIkk.lU.-i 


and 


|v“(v(v*p,)ox,)||oo<c ||v/'v(y*p,)||^||Aii“ii)| 


I ^(l"l/3|l)|| 
Ws oo 


|/3|<|a| «l,-,01/31 

oi+...+a|^l=fc 


(B.13) 


for a /c-dependent constant C > 0. Let k = 3. We have, for |a| = 3: 

||V“(V(y*p,)oX,)||oo 

|V4(y * p,)||oo||XW||^ + ||V3(y * p.)||oo||X(2)||^||xW||^ + ||V2(y * p,)||oo||Xi3)||^ 


< c 


(B.14) 


<Ce^l^l||Wo||„2+C'||Ai3)||, 
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where the constant C > 0 is allowed to depend on || j but not on the higher Sobolev norms. The 

last step follows from ()B.12p and from the previous estimates on ||bI4||„j, ||Xp^||oo, j = 1,2. Plugging 
this bound in (IB.3p . we hnd 

||‘hf^||oo<Cle^'"||Wo||H| (B.15) 

where the constant C > 0 is allowed to depend on ||lTo||j:^i, but not on the higher Sobolev norms. 
Thus, proceeding as in (jB.Sp : 

\\Wt\\Hi<Ce^^'Wo\\Hi (B.16) 

where the constant C > 0 is allowed to depend on ||1 To|Ih 2, but not on the higher Sobolev norms. 
This concludes the proof of (|B.2p for k = 3. Let A; = 4. Similarly to (|B.14D . using (IB.12P together 
with the estimates for ||VFs||j:^j, j = 1,2,3, we hnd, for |a| = 4: 


\V<^iV{V*ps)oXs 


< c 


|v5(p*p,)||oo||xW||^ + l|v4(p*p,)||^||xp)||^l|x«y 


2 

CX) 


s lloo| 


+ ||v3cp*p,)||, 

<Ce^^^\\\Wo\\m + C\\xi^'>l 


||xib||^ + ||xi2)||M + ||v2(P*p,)|U||xW||, 


where the constant C > 0 is allowed to depend on ||lTo||j:/ 2 , but not on the higher Sobolev norms. 
This implies 

||#||oo<Ce^l*l||B"o|lH3 (B.17) 

where the constant C > 0 is allowed to depend on ||lTo||//|, but not on the higher Sobolev norms. 
Then, we claim that: 

\\Wt\\Ht<Ce^\^\\\Wo\\Ht (B.18) 

where the constant C > 0 is allowed to depend on ||lTo||j|^ 2 , but not on the higher Sobolev norms. In 
fact, from (|B.5p we get and from the previous estimates on ||‘bp^||oo, j < 4, we have: 


I VP* 


1^4 < 


llbPoll 


Hi 


1^, 


(4)||2 


+ 


4 

E 

k=2 




Hi 


(B.19) 


This, together with (|B.17I) . implies (IB.lSp and concludes the proof of (|B.2p for A: = 4. The case k = 5 
can be studied in a similar way. Let |a| = 5. Using once more (IB.12p . (IB.13h . and proceeding as for 
the previous cases, we get: 


||V“(V(U * ps) o X,)||oo < C'e^l‘l||IUo||//4 + C||Xp)||^ (B.20) 

where the constant U > 0 is allowed to depend on ||IUo||j|^ 2 , but not on the higher Sobolev norms. By 

Gronwall’s lemma, we get: 

M^^\\oo<Ce^'^Wo\\Hi (B.21) 

where the constant C > 0 is allowed to depend on ||ITo 11/^2, but not on the higher Sobolev norms. 
Then, we claim that: 

||IUi||^5<Ce^l^l||IUo|lH| (B.22) 
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where the constant C > 0 is allowed to depend on HWoIIh^) but not on the higher Sobolev norms. To 
see this, we use again ()B.5p . We get: 






(5)||2 




+ 




k=3 


2k 

Hi 


$ 


(4)||2 ll^(l)||2 


0 


0 


+ 11^ 


(3)||2 ||^(2)||2 


0 


0 


(B.23) 


which, together with (jB.15p . (IB.lTp . pB.2ip . implies (|B.22p . This concludes the proof of ()B.2p for 
k = 5, and of Proposition lB.il □ 


C Propagation of commutator bounds along the Hartree dynamics 


Bounds for norms of commutators of the form [x,ujN^t] and [eV,a;Ar^t] play an important role in our 
analysis. In this section, we show how they can be propagated along the Hartree evolution. Similar 
bounds have been proven in [7]. 

Proposition C.l. Assume 

j I^(P)I(1 + \p\‘^)dp < oo (C.l) 

Let be the solution of the nonlinear Hartree equation 

iedtUJN,t = [-e^A + (H * pt),ujN,t] 

with initial data ujN,t=o = ^N- Then there exists a constant C > 0 such that 

||[x, [||[x, watJIIfs' + ||[eV,a;Ar]||FS'] 

||[eV, a;Ar,t]||F5 < [||[x, watJUfs' + ||[eV,a;Ar]||/fs'] 

Moreover, 

\\[x,ujN,t]\\tr < [||[x,a;Ar]||fr+ ||[eV, wat]|| tr] 

||[eV,cjAr,f]||tr < [||[x,cjAr]||fr.+ ||[eV,a;Ar]||tr] 

Proof. Let hnit) = —e^A + (H * Pt){x) and U{t-,s) be the unitary evolution generated by hnit), as 
defined in (13.111 . We compute 


isdtl(*{t]0)[x,ujN,tp({t;0) = -U*{t;0)[hH{t), [x,u}N,t]p({t;0) +U*{t;0)[x, [hH{t),ujN,t]p({t;0) 

= U*{t]0)[[h{t),x],ojN,tp{t-,0) 

= eU*it-,0)[eV,u:N,tp{t;0) 


Integrating over time, we find 

[x,u}]\f^t] =L({t;0)[x,ujN]ld*{t]0) + i f dsU{t-,s)[eV,uj]\f^s]h(*{t;s) 

Jo 


and thus 


||[x,a;Ar^t]||HS < ||[a;, watJIIhs + / ||[eV, wat^sJIIhs 

Jo 


(C.2) 
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On the other hand 


i£dtU*{t]^)[eV= -U*{t]{))[hH{t), [£V,u}N,tW{t;0) +U*{t;0)[£V, [hH{t),u]N,t]p{t-,0) 

= U*{t- 0) [[eV, hn (t)], ujN,tp( (i; 0) 

= £U*{t-,0)[V{V * pt),uJN,t]U{t-,0) 

Using the identity 

Jo 

we obtain, with (jC.ip . 


||[eV,a;Ar,t]||HS < ||[eV,a;Ar]||HS + J dp\V{p)\\p\‘^\pt{p)\ (is||[x, WAf,s]|| 

< ||[eV, watJIIhs + C* [ ds ||[x, WAr,s]||HS 

Jo 

Combining the last equation with ()C.2p and applying Gronwall’s lemma, we find 

[||[a;!‘^Ar,t]||HS + ||[eV,a;Ar,t]||Hs] < C'e'^1‘1 [||[x,WAr]||HS + ||[eV,a;Ar]||HS 


HS 


(C.3) 


as claimed. In the same way, one can also prove the estimates for the trace norms of the commutators. 

□ 
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